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CRYSTALLOGRAPHY 


CRYSTALLOGRAPHIC STUDIES ON SOME “COENZYMES Q” 
BY 


CAROLINE H. MacGILLAVRY 


(Communicated at the meeting of September 26, 1959) 


1. A few years ago, Dr. R. L. Lester of Wisconsin University sent 
me for X-ray investigation a number of samples of substances, isolated 
from various sources, to which the collective name of coenzymes Q has 
been given (LESTER c.s. 1958). At the time, little was known about these 
compounds, and in particular, their molecular weights were still doubtful. 

The crystals are orange-yellow coloured, of tabular habit and irregular 
shape; two prevalent edges are nearly perpendicular. 

Although the quality of the tiny crystals was mostly poor, it was 
possible to make some single crystal X-ray diagrams in a Weissenberg 
camera, from which the cell constants could be determined. The crystals 
were found to be triclinic, pseudo-monoclinic, and frequently twinned. 
Densities were measured by the flotation method and from these data 
and the cell constants the molecular weights of the substances were 
found, assuming two molecules per unit cell. Preliminary values of these 
molecular weights have been reported by LESTER e.a. (1958, 1959). The 
cell constants have since been refined by least-squares; the results are 
given in Table I. 


TABLE I 
Cell constants, densities and molecular weights of Coenzymes @ from different 
sources 
el 
Beef Heart 1 Azotobacter 
mitochondria Torolasattie vinelandii 
oe 
a 6.19 + 0.08 A 6.2 +0.15 A 6.18 + 0.08 A 
b. Be 04 A 54.85 +0.9 A 50.1 40.6 A 
c 7.8 -- O01 A 7.78 + 0.04 A 7.79 + 0.04 A 
do1o 56.9 +0.4 A 52.45 +0.9 A Alyy AEG) VAN 
loa (Ouse 91+ 1° 91 + 1° 
Be (89°) $9 41° BO. 1° 
gr emas 73 + 1° (1B) 73 + 3° 
density 1.09 1.085 1.083 gr/em? 
mojecular 
weights: 
M caps 900 820 750 
JEAN 863 795 727 


* values between brackets are by 
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analogy rather than actually measured. 


264 


(Since doi is by far the largest spacing, the b-axis could have been 
chosen in such a way that the cell would be more nearly an orthogonal 
parallelepiped. The axes in Table I were chosen so as to obtain as close 
a relation as possible to the subcell discussed below.) 

It is seen that the molecular weights thus determined are consistently 
somewhat higher than the values calculated from the structure formula I 
which has since been established for coenzyme Q (WOLF e.a. 1958, MORTON 
e.a. 1958). 


(1) | CH; 
CH;0’__) (CHe-CH=C !|— CHe),—H. 
O 
In this formula, n, the number of isoprene groups in the long chain, 
can have values varying from 6 to 10, dependent on the source from 
which the coenzyme is isolated. For the substances described in this 
study, the following values for n have been given. 


Beef heart mitochondria. . .... n= 10 
Torulasutilis .. .. ss -2«s 2. B= § 
Azotobacter vinelandii. ....s*. n= 8. 


2. From the cell constants of the three compounds studied by us, 
it is seen that their crystal structures must be very similar. This is also 
evident directly from the X-ray diagrams which show similar very 
characteristic intensity distributions. Moreover, all films indicate that 
the contents of the unit cells have a very pronounced pseudo-periodicity. 
The clusters of strong reflections define a ‘‘subcell’’ (ParrERsoN, 1931, 
VaNnD, 1951) which bears the following relation to the cells defined in 
Table I. 


subcell true cell 

{100} [100] 

[010] [021]/p, where p is such that bg poey is about 4.8 A. 
[OOL| {O01} 

(100) (100) 

(O10) (0, p/2, 0) 

(001) (0, —1/2, 1) 


This subcell, which is common to the three compounds, reveals the 
periodicity of their isoprene chains. It can thus be concluded that practi- 
cally all the isoprene units have the same configuration and even con- 
formation. 

It is interesting to compare the dimensions of this subcell with those 
of 6-guttapercha (BuNN, 1942), Table II. 

The very close relation of the cell dimensions suggests a stacking of 
the same elements. This would mean that the isoprene chains in coenzyme 
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TABLE II 


Comparison of cell dimensions ete. of 
ee 


| Coenzyme Q : 
| Subeall 6 guttapercha 
——— ee Ee See eee 
a. 6.18 A UUs) AN 
dioo 5.92 A | 5.89 A (1/2 a) 
b | 4.8 A | IN 
c | ToT) IN | 7.78 A 
Ce Set a Sie 90° 
aa he ee | 89° 90° 
oR pe Se cence. oe oe 
space group . IPT (3) | P2121 24 
number of 
molecules 
per cell 2 4 


Fig. 1. a. Structure of f guttapercha according to 
Jerrrey (1944), with origin shifted so as to conform with 6. 
b. Proposed structure for the subcell of coenzyme Q. 
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Q have the same configuration and conformation as those in f-gutta- 
percha, i.e. all-trans across the double bonds and non-planar with respect 
to the single bonds. Taking the coordinates of the atoms in the /-gutta 
unit as given by JEFFREY (1944), a tentative structure for the subcell 
of coenzyme Q was set up (fig. la and 6), and subcell structure factors 
were calculated for a number of reciprocal nodes corresponding to this 
subcell. The agreement with experiment is only fair, as could be expected 
since the whole quinone part of the molecule is neglected in these calcula- 
tions. But it seems that the model could be approximately correct for 
the “tails” of the molecules. In any case, there is little doubt about the 
all-trans configuration in the isoprene chains of these compounds. It is 
known that in rubber, with cis-configuration, the overall length of the 
isoprene unit is markedly shorter than 4.7 A. It is suggested that an X-ray 
investigation of the compounds described as ubiquinones by MORTON e.a. 
(see for example Morton, 1958a) might reveal whether they have the 
same chain-configuration as the coenzymes Q. With respect to this it may 
be of interest that a “ubiquinone” prepared by Professor E. C. SLATER 
of this University gave a powder diagram undistinguishable from that 
of the coenzyme Q from beef heart mitochondria. 


3. Haperimental details (with A. KREUGER). 


Q 275 from beef heart!). Crystals repeatedly twinned on (010). 
Abnormal interference figure under polarising microscope. Weissenberg 
diagrams of zero and first layer about axis of 7.8 A (c-axis). 

M 4 from beef heart. Less frequently twinned, unmistakably triclinic. 
Identity with Q 275 was proved by comparison of (010) reflections on 
Weissenberg diagram, and also of powder diagrams made in Guinier— 
de Wolff camera. 

Weissenberg diagrams as for Q 275. The long spacing do19 was determined 
from an oscillation diagram flanked by two Debye Scherrer diagrams of 
Al-powder for calibration. 

1-30-4 B 4 from Torula utilis. Weissenberg diagrams of only slightly 
twinned crystal about axis 6.2 A (a-axis), zero and first layer line. Optics 
the same as next sample. 

3-4-3 L from Azotobacter vinelandii. By far the best crystals, but 
still very tiny. Optics: Vibration directions in the plane (010) at an angle 
of about 17° to the a-axis and about 16° to the c-axis respectively. 
Interference figure: one binormal at the border, n, within the field of 
vision (aperture 0.85), at azimuths of 80° and —9° respectively from the 
a-axis. From these data it is possible to construct the axial plane. 
(010): n, = 33°, (010): Bn'= 37° approximately. Axial angle 2V about 90°. 
X-rays: Weissenberg diagrams about axis 6.18 A, zero and first layer line. 

The X-ray apparatus used was put at the disposal of the author by the 
Netherlands Organisation for Pure Research (Z.W.O.) 


') Sample numbers as given by R. L. Luster. 
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PHYSICAL CHEMISTRY 


INFLUENCE OF THE IRRADIATION TEMPERATURE ON THE 
POSITION OF THE F-BAND IN MIX-CRYSTALS OF KCl AND NaCl 


BY 


BE. OTTENS, A. J. ELAND!) anny W. G. BURGERS 
(Communicated at the meeting of September 26, 1959) 


Summary 

The position of the F-band is measured for mix-crystals of NaCl and 
KCl, obtained by quenching from the homogeneous state. Measurements 
are made on crystals X-rayed both at room temperature and at the 
temperature of liquid air. Also the absorption is measured at both temper- 
atures. A special attachment for the Beckman-spectrophotometer is 
described, which enables such measurements to be carried out. 

It was found that for the mix-crystals the position of the F-band is 
dependent on the temperature at which the crystals are irradiated, the 
maximum lying at shorter wave-length (difference about LO-15 my) for 
irradiation at liquid air temperature as compared to irradiation at room 
temperature. Such a shift is not observed for the pure crystals: here the 
position of the band is independent of the irradiation temperature. 

For mix-crystals irradiated at liquid air temperature the F-band lies 
between that for pure NaCl and pure KCl. The ‘‘anomalous’’ position of 
the maximum of a 10 NaCl/90 KCl crystal irradiated at room temperature 
as observed by Gnaedinger (position at a longer wave-length than that 
for pure KCl) disappears for irradiation at liquid air temperature. 


1. Introduction 


For some years the decomposition of mix-erystals of NaCl and KCl 
has been studied in this laboratory. It is well-known that below 500° C 
the mutual solubility is limited and decreases with decreasing temperature. 
X-ray and thermodynamic investigations concerning the phenomena 
accompanying the decomposition were carried out by several investigators, 
in particular by TicHELAAR (1956), who also discussed literature up to 
that time. 

We have started an examination of F-centers in mix-crystals, because 
such centers may be expected to be sensitive to changes in ion con- 
figuration in their immediate neighbourhood 2), On the other hand it 


') Now at Philips Research Laboratories, Eindhoven. 
*) See e.g. the theoretical considerations of GNAEDINGER (1953) on the “structure” 
of F-bands, in which the author considered them to be built up of seven “elementary” 
bands, due to halogen vacancies with different surroundings, viz. 6 Na*, 5 Nat 4-1 Ke, 
4 Nae + 2Ks,...6K-, 
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may be concluded from Tichelaar’s work that at temperatures slightly 
under the solubility line the decomposition takes place by displacements 
of the cations without essentially breaking up the single crystal structure. 
Tn its first stages such a process would involve an enrichment in number 
of equal ions (either Na or K) in small lattice regions !). This would 
influence the energy levels of the electrons in the F-centers and therefore 
the position of the corresponding F-bands. For this reason it seems 
possible that the behaviour of the F-bands might lead to a closer under- 
standing of the decomposition mechanism. 

Up to now, however, we have only studied the formation and the 
position of F-bands in the homogeneous mix-crystals. The only author, 
as far as we know, who reported on the subject, is GNAEDINGER (1953). 
In this paper measurements are given concerning the position of the 
F-band in mix-crystals of KCl and RbCl, which salts form a continuous 
series of stable solid solutions. It was found that the positions of the 
maxima were intermediate between those of the pure components. This 
is what one would expect on account of a rule given by MoLiwo (1933), 
according to which the position of the maximum of the F-band, at a 
given temperature of the crystal, only depends on the lattice constant d 
according to the relation /,,,, =const. d?. Gnaedinger’s paper also contains 
measurements concerning a metastable quenched mix-crystal of NaCl 
and KCl of composition 10 NaCl/90 KCl (mol. %). Most remarkably for 
this crystal, the F-band was situated on the “wrong” side of the F-band 
of pure KCl, viz. at a longer wave-length and therefore outside the range 
between the F-band of pure NaCl and pure KCI. 

It should be noted that in these experiments the crystals were X-rayed 
at room temperature. We thought it of interest to irradiate also at lower 
temperatures. For this purpose a special attachment to a Beckman- 
spectrophotometer was constructed, the description of which is given in 
the following section. 


2. Haxperimental 
The spectrophotometer used was a Beckman-instrument, model DU. 
The cell-holder of this instrument was replaced by a special attachment 
(cf. fig. 1) in order to adapt the instrument to our measurements. This 
apparatus had to meet the following requirements: 
a) Irradiation and measurement had to be possible at any temperature 
between that of liquid air and of about 150° C. 
b) Condensation of water vapour on the crystal had to be avoided 
absolutely, because even traces of water have a decomposing effect 
on the metastable mix-crystals. 


1) Indications based on X-ray diffraction effects, obtained by KioorwisK and. 
TrrpEMA (1958), made it probable that a certain degree of ordening is present in 
the distribution of the Na- and K-ions in a decomposed 50-50 crystal. 


This necessitated irradiation of the crystal in vacuum. In fig. 2 a 
sectional view of the apparatus is given. The vacuum cryostat consists 
of a stainless steel dewar flask A, mounted on a square metal plate B. 
By means of this plate it is possible to position the cryostat accurately 
in the photometer as well as in the irradiation-house attached to the 
X-ray tube. 

The original cell compartment of the photometer is replaced by the 
metal box C. As will be described in section 3, for each measuring-point 
it is necessary to adjust the intensity of the incident light-beam. For this 
purpose the sample has to be taken out of the light-beam and replaced 
by a calibrating opening. To make this possible the cryostat is fastened to 
a carriage D, rolling on a number of steel balls E lying in “V” grooves. 
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E. OTTENS, A. J. ELAND anp W. G. BURGERS: Influence of the irradiation 


temperature on the position of the F-band in mix-crystals of KOL and NaCl 


Fig. 1. Attachment adapted to the Beckman-spectrophotometer for X-ray 
irradiation and absorption measurement of crystals at various temperatures. 
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The slide is pressed down on the steel balls by two spiral springs F. In 
this way a “tumbler-action” is obtained giving two reproducible positions 
bringing either the crystal or the calibrating opening in the light-beam. 
At the bottom the dewar flask is equipped with two quartz windows G, 
allowing light transmission through the crystal. The holder with the 
crystal (cf. fig. 3) is inserted in opening I of the copper house K and 
locked with a small screw. 


Fig. 3. Crystal-holder. 


The crystal can be irradiated through the plastic-foil window H by 
rotating the crystal-holder 90 degrees by means of the joint at L (sealed 
by an “O” ring). 

The copper house K is fixed to the bottom of the liquid-air container 
by means of a stainless steel tube M, while, in a straight line with M, a 
brass tube N is fitted. By inserting a copper rod in this tube down to the 
bottom, a good heat-exchange is obtained between the copper house K 
and the liquid air. By replacing the copper rod by a small electric furnace 
(a glass tube containing a small heating-element at the bottom) any 
temperature from —180°C to about 150°C can be obtained in a few 
minutes. The low thermal conductivity of the stainless steel tube M 
keeps heat-loss down. E.g. only 5 watts are required for maintaining a 
temperature of 0° C. 

The temperature was measured by a copper-constantane thermo- 
couple, the leads of which were brought out uninterrupted through 
kovar seals with hollow inlets, thus preventing undesired thermal poten- 
tials. Because the photometer was fitted out with a photomultiplyer 
tube, which was light-sensitive to a high degree, it was important to keep 
out any light. This was achieved by enveloping slide D, as it were, in a 
flat box O having a play of a few tenths of mms. Moreover all internal 
parts were painted black with camera lacquer. 

The inner mantle of the dewar flask was made as thin as possible at 
its top in order to keep heat-supply through the metal cylinder as low as 
possible. Nevertheless a circular chamber P, containing a small heating- 
element, was necessary to keep the temperature of the top of the dewar 


flask above 0° C. 
The crystal-holder (cf. fig. 3) consists of two brass blocks A and B, each 
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containing two openings, 1 and 2, of equal diameters. The crystal plate 
is laid upon opening 2. Then block A is placed as a cover on B. The crystal 
is fixed by two small springs C. In order to insert the crystal-holder in 
the copper house K (ef. fig. 2), the dewar is dismounted at L. 


3. Measurement 


As we wanted to measure the absorption caused by the F-centers, the 
absorption of the non-irradiated crystal had to be subtracted from the 
absorption of the coloured crystal. To achieve this for a chosen wave- 
length, slide D is first placed in the position in which the light passes 
through opening | of the specimen holder (intensity Ip = 100). D is then 
shifted to its other position in which the light passes opening 2, and thus 
passes through the crystal (intensity [). In this way the absorption, given 
by the optical density D-=In Io/I, is determined both for the non- 
irradiated and for the irradiated crystal. At the same time a possible 
difference in the diameters of the openings | and 2 is eliminated from 
the measurement. 

The X-radiation originated from a tungsten tube, run at 47 KV and 
17mA. The duration of each exposure was 30 minutes. The crystal 
plates had dimensions. of approximately 4 x 3 0.5 mm, 


4. Preparation of crystals 


The mix-crystals were prepared in the same way as described by 
TICHELAAR (1956). Known quantities of NaCl and KCl (Brocades p.a.) 
were mixed in the pulverized state and then melted in a porcelain crucible 
in an electric furnace. The melt was cooled slowly to a temperature close 
below the solidus-line and kept there for about two days to homogenize. 
It was then cooled slowly to about 550° C (the top of the solubility curve 
lies at about 500°C) and finally quenched to room temperature. The 
cleaved parts thus obtained were single crystals; they were kept in an 
exsiccator with P2O5. All further manipulations were carried out in dry 
air, to prevent decomposition. For a short time the crystals can be pro- 
tected against humidity by immersion in zapon varnish, diluted with 
amylacetate, by which process they become covered by a thin layer of 
lacquer. The compositions were determined by X-rays, using the relation- 
ship between lattice parameter and composition as given by Bunk and 
TICHELAAR (1953). Crystals were prepared with the following compositions 
(in mol. % NaCl-KCl): 100/0, 88/12, 65/35, 53/47, 39/61, 18/82, 9/91, 0/100. 


5. Results 


Distinct bands, lying in the range of wave-lengths expected for F- 
bands, were observed for mix-crystals with compositions not too far from 
the pure components, namely for 9/91, 18/82 and 88/12 NaClKCl. 
Crystals with compositions nearer to the aequimolecular composition 


Orn 
Dts 


either did not show such bands or they were very unstable !) 2). 
. As an example we show in fig. 4 the absorption spectrum for the compo- 
sition 88 NaCl/12 KCl. The figure shows four curves, obtained with two 
crystals. Crystal I was X-rayed at room temperature (about 20° C) 
crystal IT at —180°C. For both crystals the absorption was nonce 
at 180" C and at.20° C. 

The figure first of all shows the well-known fact, characteristic for an 
F-band, that its maximum shifts to a longer wave-length with increasing 


) 


To 
I 


—— Optical density (log 


—— Wavelength in my 


Fig. 4. Absorption spectrum of two X-rayed mix-crystals with composition 
88 NaCl/12 KCl. Crystal I: X-rayed at 20°C. Crystal Il: X-rayed at —180° C. 
The absorption of both crystals was measured at ——-180°C and at 20°C. 


1) For the system KCl-RbCl it was found by GNAEDINGER (1953) that the more 
the composition approaches 50-50, the larger the half-width of the F-band. 

2) Moreover, most crystals showed an “abnormal” band at about 420 my, 
discussed by Hanperin and ScHLESINGER (1959) and by OrrEns et al. (1959). 
As shown in the last-mentioned paper, the intensity of this band was lower when 


irradiated at room temperature. 
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temperature of the absorption measurement, while at the same time 
the band-width increases. The first fact is a direct consequence of 
Mollwo’s rule (cf. section 1) in connection with the thermal dilation of 
the crystal lattice, the second is due to the increased thermal lattice 
vibrations. 

Moreover it appears that the position of the maximum is dependent 
on the temperature at which the crystal is irradiated, the maximum 
lying at a shorter wave-length for the crystal X-rayed at —180°C as 
compared with the crystal X-rayed at room temperature. The shift due 
to this cause amounts to about 10-15 my. A similar dependence of the 
maximum on the temperature of irradiation is not found for the crystals 
of the pure components and it is therefore characteristic for the solid 
solution state. 

In fig. 5 a graph is given of all the maxima measured as a function of 
composition, and the data are put together as well in Table I. 


TABLE I 


Wave-length of the maxima in the F-bands in NaCl-KCl mix-crystals as function 
of composition, temperature of irradiation and temperature of measurement 


, Wg: Measured at 
siete 3 da es X-rayed at 
mol. % NaCl 20° —180°C 
100 (NaCl) 20° C i a 
180° C 462 450 
88 20° C 495 480 
—180° C 485 470 
18 20° C 565 540 
—180° C 550 532 
9 20°. C 567 544 
—180° C 552 534 
0 (KCl) 20° C . 
180°C 557 540 


6. Discussion 


An inspection of fig. 5 and Table I learns that for crystals X-rayed at 
liquid air temperature the positions of the maxima found for the mix- 
crystals lie between those of the pure components. This holds good in- 
dependent of the temperature of measurement. The “anomaly” observed 
by Gnaedinger (ef. section 1) does not occur for irradiation at low tempera- 
ture. Our results confirm Gnaedinger’s finding for crystals irradiated at 
room temperature: under these circumstances both the crystals with 
composition 9 NaCl/91 KCl and those with 18 NaCl/82 KCl have an F-band 
with a maximum wave-length larger than that belonging to pure KCl. 

It seems therefore that irradiation at liquid air temperature produces 
a different state in a mix-crystal, possibly different surroundings of the 


Measured at 20°C 
A X-rayed at 20°C 
590)—=\ == = as 160°C 
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Mollwo's rule 
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Fig. 5. Position of the maxima in the F-bands in X-rayed NaCl-KCl mix-crystals as 
function of composition, temperature of irradiation and temperature of measurements. 


F-center, as compared with the state produced by irradiation at higher 
temperature !). As mentioned in section 1, it might be expected that 
such a change influences the position of the F-band. 

As already stated in section 1, an intermediate position of the F-band 
maxima for the mix-crystals, relative to that of the components, can be 
expected on account of Mollwo’s rule in connection with the continuous 
increase of lattice constant from pure NaCl to pure KCI. If this rule would 
be the only determining factor, then, considering the linear increase of 
the lattice constant and the fact that the total change amounts to only 
10 %, one might expect that in first approximation the position of the 
maximum of the F-band would also vary in a linear way with the compo- 
sition (ef. fig. 5). It is evident that the actually observed maxima do not 
coincide with the calculated ones: the former all lie at longer wave- 
lengths. This again is an indication that for the mix-crystals the lattice 
constant is not the only determining factor and points anew to an influence 
of the actual surroundings of the F-centers. 

In this connection it is of interest to notice that Mollwo’s rule cor- 


1) In a discussion with Professor W. DEKEySER of Ghent this possibility was 
raised by him. 
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responds fairly well with the change in position of the F-maximum found 
for various pure alkali salts at normal concentrations of F-centers. 
However, measurements of KatsEr (1952) on evaporated salt films show 
that for excessive concentrations of halogen-vacancies the F-band shifts 
over about 120A (12 my) in the direction of shorter wave-length. 
According to Srrtz (1954) one may assume in such a case that already 
in the first shell of halogen-ions, surrounding a given halogen-vacancy, a 
halogen-ion is missing. The surroundings of each F-center would be dis- 
turbed in this way and this might be the cause of the observed shift. 

From fig. 4 it may be concluded that the percentage of F-centers, 
disappearing by thermal bleaching (by increasing the temperature of the 
crystal from liquid air temperature up to room temperature), is much 
larger for the crystal X-rayed at — 180° C (crystal IT) than for the crystal 
X-rayed at 20° C (crystal I). As it may be assumed that the disappearance 
is mainly due to the recombination with Vi-centers (positive ion-vacancy 
plus an electron-hole) and as the number of recombinations will be 
proportional to the number of Vj-centers present, it seems probable that 
irradiation at — 180° C produces more V,-centers than irradiation at room 
temperature. This is actually observed for KCl crystals (DORENDORF 
(1950) (1951)). For mix-erystals the “shoulder” (a) in fig. 4 at about 
360 mu points to the same conclusion. 


Laboratory of Physical Chemistry, 
Technological University, Delft, 
The Netherlands 
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ABSORPTION BAND AT 420 mu IN X-RAY COLOURED 
KCLCRYSTALS AND IN MIX-CRYSTALS OF KCl AND NaCl 


BY 


EK. OTTENS, A. J. ELAND 1), H. B. ZEEDIJK anp W. G. BURGERS 


(Communicated at the meeting of September 26, 1959) 


Summary 

In X-rayed crystals of KCl and in mix-crystals of NaCl and KCl an 
absorption band has been observed at a wave-length of about 420 my. 
The presence of this band is probably due to some unknown impurity 
in the crystals. The measurements confirm the observation of HALPERIN 
and SCHLESINGER (1959) that the intensity of the band is highly dependent 
on the thermal pre-treatment of the crystals; moreover it is unstable at 
room temperature. 


1. Recently HALPERIN and SCHLESINGER (1959) reported the appear- 
ance of a thus far not-observed absorption band at about 420 my in 
X-rayed KCl-crystals which had been heat-treated (kept at 500° C during 
several hours) before irradiation. The intensity of the band, as compared 
with the normal F-band, was markedly dependent on the thermal pre- 
treatment. The band could be bleached by irradiation at its own wave- 
length under formation of additional F-centers, and vice versa. Thermal 
bleaching has been observed as well. According to the authors, a correlation 
seems to exist between the new absorption band and the thermolumi- 
nescence of the crystals. 


2. Apparently we had come across the same band in the course of an 
investigation of F-centers in X-ray coloured mix-crystals of NaCl and 
KCl2). As an example fig. 1 shows the absorption of a mix-crystal with 
composition (in mol. %) 9NaCl/91KCl, obtained after melting the compo- 
nents together in a porcelain crucible and after quenching from a tempera- 
ture above the solubility curve (cf. TicHELAAR (1956). The crystal was 
irradiated and the absorption measured at — 180° C. In this case a maxi- 
mum is present at about 420 my, which is much higher than that at 
530 my, due to the F-centers. 


1) Now at Philips Research Laboratory at Eindhoven. 
2) Of. Orrens, ELanp and Buresrs (1959). In this paper also details of the 


experimental method are given. 


100 


9NaCl-91KCI 
X-rayed at —180°C 


Measured at - 180°C 


° 

al Loe 

o 080 
2 

— 

3 

c 

® 

U 

8 060 
S 

| 040 


020 


0 
300 400 500 600 


——e Wavelength in mE 


Fig. 1. Absorption-curve of mix-crystal 9 NaCl/91 KCl (mol. %). X-rayed and 
measured at 180° C. Strong 420 my band. 


As mix-crystals with various other compositions, prepared from the 
same material (Brocades p.a.), showed a band at approximately the same 
wave-length, we thought at first that its occurrence was characteristic 
for the solid solution state. This, however, appeared not to be the case. 
The band was also observed in crystals of KCl used for the preparation 
of the mix-crystals. This is shown in fig. 2, again for irradiation and 
measurement at —180°C, The two curves correspond to a) a crystal 
quenched from 550°C and b) a slowly cooled crystal. The pronounced 
difference between the two curves confirms the conclusions of Halperin 
and Schlesinger that the intensity of the band is strongly dependent on 
the thermal pre-treatment of the crystals. It also confirms the existence 
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Fig. 2. Absorption-curve of two X-rayed KCl-crystals. X-rayed and measured 
at —180° C. Curve a: erystal quenched from 550° C: 420 my band strong, F-band 
weak. Curve 6: crystal slowly cooled: 420 my band weak, F-band strong. 


of an interrelation between the 420 my-band and the F-band, in the 
sense that the one is stronger in case the other is weaker. 


3. In the course of our work we had also obtained KCl-crystals 
(prepared from material of unknown origin) which did not show a trace 
of the new band, although the crystals were quenched from 500° C. This 
raised the question whether its presence was solely an effect of thermal 
pre-treatment, or whether also the purity of the material played a part. 
In our method of preparation the possibility could not be wholly excluded 
that some ‘impurity’? might have been introduced into the melt. We 
therefore prepared crystals (from KCl, Brocades p.a.) under such pre- 
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Fig. 3. Absorption-curve of two “‘pure’’ KCl-erystals cleared out of the same 
block and both heated at 500° C for 5 hours. X-rayed and measured at —180° C. 
@: crystal slowly cooled. ©: erystal quenched. 420 my band absent. 


cautions as would prevent or at least reduce to a minimum the occurrence 
of an impurity, e.g. the porcelain crucible in which the salt was melted 
was cleaned beforehand in boiling aqua regia, the crucible was covered 
during melting and the cover was not once lifted in the course of the 
preparation etc. Crystals prepared in this way did not show the 420 my- 
band, neither could it be produced afterwards by a heat-treatment, as 
apphed by HaLprrin and ScuiesincEer. This is obvious from fig. 3, 
which shows the absorption curves of two crystals cleaved out of the 
same block. Both crystals were heated at 500°C for 5 hours, then one 
(@) slowly cooled, the other (0) quenched. This result strengthens our 
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Fig. 4. Influence of temperature of irradiation on the 420 my band of a KCI- 

crystal. Curve a: X-rayed and measured at —-180° C. Crystal then bleached at 

100° C for 4 hour. Curve 6: X-rayed at 20° C, measured at —180° C. 420 my band 
weaker, F-band stronger. 


supposition that the band is primarily due to some “impurity” in the 
KCl used by us. At present we do not know its nature. In a recent paper 
Isaiguro and OKUNO (1959) report a band at 425 my for silver-doped KCl. 
They call this the E-band. According to the authors it is due to the 
presence of neutral silver atoms. It is remarkable that also the K-band 
shows a “‘strong-weak”’ relation with regard to the F-band, as may be 
seen from Ishiguro and Okuno’s paper. 


4, Finally we refer to the observation of HALPERIN and SCHLESINGER 
that the 420 my-band can be bleached at about 250° K. This makes it 
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probable that the band is unstable at room-temperature. In fact, keeping 
an irradiated erystal at room-temperature for some time makes it disappear 
completely. Moreover, irradiation at this temperature causes a much 
weaker intensity than irradiation at liquid-air temperature, in favour of 
the F-band. This is shown in fig. 4. 


We intend to continue this investigation. 


Laboratory of Physical Chemistry, 
Technological University, Delft, 
The Netherlands 
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GEOLOGY 


AN OCCURRENCE OF CORUNDUM IN THE PRECAMBRIAN 
OF GHANA 


BY 


A. F. J. SMIT anp J. E. J. M. VAN LANDEWIJK 


(Communicated by Prof. G. H. R. von KoEnrcswa tp at the meeting of Sept. 26, 1959) 


Corundum has occasionally been found in placer deposits in Ghana 
but as far as we know there are no records of the mineral in situ in the 
country. This preliminary note records the discovery of an occurrence of 
corundum in situ, together with a brief account of the geology and 
mineralogy of the occurrence. A more detailed account will be published 
later. 

The rock containing corundum was found near Nebiewale in N. Ghana. 
The geographical position of this village is 1°40’ W. Long. 10°35’ 20” N. Lat. 
(sheet Navrongo 1944, Topographical map of Ghana 1 : 250,000) 1). 

The Nebiewale region is underlain by granite-gneisses of the so-called 
Cape Coast complex (Geological map of Ghana, 1955) which is intrusive 
into Middle Precambrian Birrimian regional metamorphic rocks, among 
others sedimentary and volcanic schists and metamorphosed basic 
intrusions. 

In the Nebiewale area granite-gneiss, cut by numerous pegmatite veins, 
predominates. The corundum occurs in lenses of muscovite-gneiss of 
varying size in the country rock. All transitions are found from clearly 
defined inclusions to ghost inclusions, only recognizable by the mineral 
composition. The rock around the inclusions is probably endomorphically 
contaminated. 

The country rock has the following mean composition: quartz 30-35 % ; 
microcline, orthoclase, perthite, 10-30 °%; albite, 0-32 %%; oligoclase, 
0-50 %; biotite, 2.5-3.5 %; muscovite, 1-1.5 %; and epidote 0.5-1 %. 
If the country rock can be regarded as an igneous rock, it is a leuco- 
sodaclase adamellite (17” P), according to JoHANSSEN. Most of the rock 
is porphyritic and heterogeneous, but the composition varies normally 
within the above-mentioned limits. It is possible that the enclosing rock 
is a migmatite, probably in contact with an intrusion in the N.E. 


1) In the literature Nebiewale has been mentioned in a paper by ROHLEDER 
who believes that the village is situated in an old volcanic crater, which he called 


the Nebiewale Caldera. 
We intend to discuss this problem in a separate publication in due course. 
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The quartz is partly myrmekitic, microcline as well as isomicrocline is 
probably present, the plagioclase is untwinned or albite-twinned. The 
biotite is pleochroic light to dark brown with many pleochroic haloes, 
some with a zircon core. The myrmekite seems to replace the plagioclase, 
the microcline the quartz; the muscovite seems to replace the biotite and 
the plagioclase. Dark inclusions, sometimes with corundum, and veins of 
pegmatite and aplite are present. 

Part of the enclosing rock has a quartz content of 10-20 %, and an 
albite content of 50-60 °%. This rock is probably contaminated by 
digestion of the inclusions. 

In parts where corundum remnants are still found, the average compo- 
sition is: quartz, 35-40 %; oligoclase, 40 °%,; biotite 20-25 %,. 

The endomorphose seems to be: biotite, oligoclase and quartz giving 
rise to the formation of alkali feldspars and muscovite, which evolve to 
quartz, alkali feldspars and oligoclase. 

The nearly-digested corundum inclusions contain andesine and a high 
percentage of biotite and no alkali feldspars. The corundum is partly 
altered into small muscovite flakes, lying at random. which form larger 
crystals at some distance from the corundum. 

Some corundum crystals are poikilitic, enclosing small specks of biotite 
oriented in one direction. The biotite seems to be older than the corundum 
which, in its turn, is altered to muscovite. 

The biotite around the corundum has no pleochroic haloes and is partly 
smaller and not so well oriented as in the rest of the inclusions. The 
clearly bounded corundum-bearing inclusions, up to 6 sq. feet, have a 
white zone of one foot around them, mainly consisting of feldspars (quartz, 
0.5 %; oligoclase-andesine, 91 °; biotite, 3.5 °; muscovite, 5 %, some 
corundum relies). 

A zone rich in biotite or composed of nearly-digested biotite-rich 
inclusions sometimes surrounds the white feldspathic rim. 

The inclusion itself is greenish in colour and has abundant, large crystals 
of grey, speckled, barrel-shaped emery and opaque to translucent corun- 
dum crystals (up to 2”), which are fluorescent, violet, red-violet and 
grey-blue. 

The corundum is generally arranged in lines, parallel to the orientation 
of the inclusions (biotite-muscovite). The crystals sometimes contain 
inclusions of iron ore or biotite and are partly replaced by muscovite as 
outlined above. 

The corundum is sometimes uniaxial, sometimes biaxial; it has a good 
cleavage and shows polygonal growth lines around a more massive core. 
The mineral is sometimes twinned. 

Muscovite, iron ore, oligoclase-andesine and some small biotite crystals 
accompany the corundum. 

We acknowledge the financial support from the University College of 
Ghana which enables us to carry out this research, 
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MECHANICS 


APPLICATION OF LIFTING SURFACE THEORY TO SHIP SCREWS 
BY 


J. A. SPARENBERG 


(Communicated by Prof. W. P. A. vAN LAMMEREN at the meeting of Sept. 26, 1959) 


Summary 

The problem of determining the relation between the geometrical form 
of a propeller blade and its pressure distribution is reduced to a two 
dimensional integral equation. This has been done by extending the 
existing lifting surface theory for airplane wings to rotating helicoidal 
surfaces. The essential tool is the acceleration potential. 


1. Introduction 

The theory of lifting surfaces has been applied extensively to airplane 
wings of finite span (for instance [1]!), [2] and [3]). In this article we 
discuss the adaptation of this theory to rotating screw propellers. Such a 
theory is needed especially for screws of fast ships, because the aspect 
ratio of their blades becomes rather small. 

These screws are generally designed by using the vortex line theory. 
This theory takes care of a mean value of the induced velocity at the 
profile. For a blade with a small aspect ratio however the variation of the 
values of the induced velocities along the chord are also important. In 
order to allow for this, several types of corrections (so called vortex line- 
vortex sheet corrections) are suggested ({4] and [5]). Recent measurements 
of vAN MANEN and CrowLEy [6] gave rise to some doubt about the 
prevailing types of corrections. In order to obtain an insight into the 
accuracy which is arrived at by the calculations just mentioned, numerical 
results from an exact lifting surface theory are necessary. 

Before doing this it is desirable to have at our disposal a more general 
treatment of the screw. To some extent this has already been given by 
GuILLOTON [7] and SrrscHELETzKy [8]. These authors replace the 
propeller blade by a layer of bound vortices. The lifting forces are caleu- 
lated by the law of Kutta—Joukowski. 

The point of departure used in this article is the theory of the acceler- 
ation potential as described by PRaNnptL [9]. Theoretically this yields 
a more satisfactory approach because then we start from the basic 
equations of hydrodynamics and do not need laws which proceed from 


') Numbers between square brackets refer to the bibliography at the end of 
this article, 
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these equations by integration. The given theory results into two equi- 
valent integral equations both of which can be chosen as a basis for 
further numerical calculations. 


2. Formulation of the problem 


The propeller blades are assumed to be infinitely thin and the influences 
of the ship and the propeller boss are neglected. First we consider the 
theory of a one bladed screw, which already contains all the specific 
mathematical complications. In the last paragraph of this paper reference 
is made to screws with more blades. 

The one bladed screw rotates with an angular velocity w around the 
« axis (fig. 1). The undisturbed velocity of the non-viscous and incom- 
pressible fluid has the magnitude U and is in the direction of the positive 
% axis. 


Fig. 1. Scheme for a one bladed screw. 


We consider a rotating helicoidal surface 
(2.1) H(z,7,9,t) = y+ot—az=0, a=F 


where ¢ represents the time and z, r and » form a system of cylindrical 
coordinates. The domain of x and r can be described by 


(2.2) CTEM): Treacy, 


the functions 2;(7) and x;(t) determine respectively the leading edge and 
the trailing edge of the blade, 7; and 7o are the inner and the outer radius 
of the blade. For a=a/U (2.1) the rotating surface H—0 does not disturb 
the incoming uniform flow. . 

Next we consider a screw blade #'=0 which lies in the neighbourhood 
of 7 =U, 
(2.3) P(x, 1, 9,4) x y+ot—ax-+ ef(x, r) = 0. 


The function f(z, 7), which is a measure for the deviation of the screw 
blade from the helicoidal surface, is some reasonable smooth function, 
ce is a small constant. This deviation (angle of attack and camber) induces 
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disturbances of the uniform flow, which are supposed to be small with 
respect to U. The components of the velocity of the fluid are denoted by 


(2.4) U+4u, v, w 


respectively in the x, r and direction. 
For an arbitrary surface G(x, r, p,t)=0 moving in a non-viscous and 

incompressible fluid the formula 

0G oG dG 0G 

bel — i —=0 
(2.5) ot 7 (U+u) ox % or is rop : 
holding on G=0 provides, the condition for the velocity to be tangential 
to the surface. The disturbance velocities wu, v and w have to be taken all 
at the same side of the surface. Substituting (2.3) into (2.5) and taking 
into account only first order terms in uw, v and w we obtain 


(2.6) 60 aude. 
Mv a 


This equation can also be derived by considering the relative velocity 
of the fluid with respect to the screw blade. 

Thus we need to determine velocities u, v and w which meet the equations 
of motion (discussed in the next paragraph) and equation (2.6). Besides 
this, wu, v and w have to vanish for 7 > —oo. 


3. The screw blade as a pressure dipole layer 


The equations of motion of a linearised time dependent flow in eylin- 
drical coordinates are 


: Se OY LC 2) 
vest) U we 8p ow , 

ow w 1 op(a, r, g, t) 

2 pa le ee) ee, 
ee) U Gina as tere 
ke ow dw 1 dp(a, r, g, t) 
3. badly I fate = etch Ay Ad 
ee) H ox § Ot ia rop ; 


where yw 1s the specific density and p the pressure. The equation of conti- 
nuity and the potential equation are 


(3.4) OU 1 orv Low _ 


‘ oe LPO ae 1 0? 
(3.5) (satrap tries) MON) =0. 


Adding the derivatives of r-(3.1), r-(3.2) and (3.3) with respect to a, r 
and » respectively and observing the equation of continuity, we find that 
p(x, r, p, t) satisfies the potential equation (3.5). Hence the pressure can 
be represented by a distribution of pressure poles and pressure dipoles on 
the screw blade. The dipoles have their axes normal to the surface of 
the blade. We shall show that only the pressure dipole layer is needed, 
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We consider the screw blade F =0 which slightly disturbs the incoming 
homogeneous flow with undisturbed velocity U. The unit normal to this 
surface has the components 


(3.6) (=) (32 2 Ly \"2(0F OF IF) at y (OF WF IF 
\\dz or rp) |} de’ dr’ ro (de? Wr? rd}? 


respectively in the x, r and direction. 
The normal acceleration 4, of a fluid particle in the neighbourhood 
of F=0 becomes 


( Of /;,0UuU . du of /,,v , w Of ow dw 
A =O (US 4. S45 -(U T+ )+ |= 
7) \ fi a Asa Piae: ou Ot "Ot rp Ose 7 ot 

i < 


oO OF UF OF o> Od 
a ( _—-> — +7) — wp — —_ 
/ (U ox =) ( ot : or uM =) (us et Os: a aor =) (U on i =) Ff. 


From (2.5) it follows that within terms of the first order in wu, v and w 


(3.8) As —0(U 4 =) F. 


This formula is independent of the side of the surface, it only depends 
on the point which we consider. Hence the normal acceleration A, is 
equal at both sides of the blade. From (3.1), (3.2) and (3.3) it follows 
that the same holds for the normal derivative of the pressure. This implies 
that the pressure in the fluid can be represented by a pressure dipole 
layer on the screw blade [10]. Because our theory is linearised this dipole 
layer can be taken to coincide with the helicoidal surface. 


4. The velocity field induced by a rotating pressure dipole 
A pressure pole of strength S, situated at the point £, 0, 6 has the 
pressure field 
def S def S 


(4.1) P(z,7,9,$,0,8) = 4n{(x — &)2 +1? +0? —2rp cos (6 —g)}'/2 a 4nR(x, 7,0, €, 0,6) 
Differentiation of this expression with respect to €, 9 and @ in the direction 
of the unit normal n of H =0 (fig. 1), yields the pressure field of a pressure 


dipole with axis normal to the helicoidal surface. The unit normal n has, 
in accordance with (3.6), the components. 


(4.2) n= (#45) (2,0, -=). 


Herwith the pressure field of a pressure dipole attached to the rotating 
surface H=0 becomes 


def So 0 1D = 
Ne Bese! ) See (# SE p36) POET PrE8 A)lomas—on 
ie) / = S {ao(a —&) +r sin (a& — wt — ¢)} 
in(l-argyh Rw, 7, 9, 6 0,48 —ot) 
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The external force exerted by this pressure field can be determined by 
calculating the resultant pressure on a small closed surface which surrounds 
the point of singularity. It turns out that this force has the magnitude S 
and is directed along n. 

The equations of motion (3.1), (3.2) and (3.3) determine the velocities 
which belong to the rotating pressure field (4.3). The right hand sides of 
these equations become 


(4.4) —1__ , | eee ee) 
_ woe Amp(1+ a2?) R(x, r, p, €, 0, a€ — wt) Rs |’ 
\ — L oD _ S [: 3{ao(a —£) +r sin (aé —wt —¢)}{r —e cos (af —at —9@)} 4 
wor 4atu(1 + a?9?)"/s Rs 
(4.5) 
] sin (a& — wt — ~) 
oe a 
\ 1 oD S [ 3 {ao(a — &)+rsin (a& — wt — ¢)} o sin (af — wt —¢) 
(4.6) ) erp Sxm(l a0 Be i 


(ag — wt — ¢) 
/ Pi 7 P|. 


The solution of 
(4.7) U 4 Y= Lia,r, 9, #), 
where L(x, r, py, t) is some prescribed function, which tends to zero for 


x —> —oo, is 


e—€ 


(4.8) y=% i" L(t +87, 9. t+ =F") de. 


From this equation we find for the velocity field induced by the rotating 
pressure dipole at time t=0 


" 5 o—€ 
(u, Vv, w) ae 4nuU (1 -+- a2)" | 
(4.9) Ee +r sin (a(a—t) —¢)}{r, r— 9 cos (a(x —t) ——), —@ sin (a(x —t)—q)} 
| (272+? — Bro cos (ale — 1) — gs BT 
_ {ag, sin (a(a—t) — ¢), — cos ae a hes Qs) 
\ {72+ r2+ 02 — 2re cos (a(x g)} ila | id 
This can be written as follows, 
(u,v, w) = 
S \ {—r sin (a& — 9), (w — &) sin (a& — ¢), — (w—€) cos (aE—g)} 
fapiras aro?) ( {(a — &)? + 72 + 9? — 2rq cos (af — ¢)}" > 


4, 
(4.10) ¢ ip [{are cos (a(a—t) —~) —ae2, —apt cos (a(x—t)—q) + 


—e sin (a(x—t)—¢), —r+¢ cos (a(a—t)——) + 
~ agt sin (a(@— rt) — 9) }{t? + 12+ 2 — 279 cos (a(x—t) —y)}-"*] dr 3 
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The proof of this statement is rather lengthy and needs only elementary 
analysis, therefore it will not be discussed here. 

In expression (4.9) the highest exponent of the two numerators, under 
the sign of integration is 5/2, while in expression (4.10) it is 3/2. This 
mathematical reduction has also physical meaning. From the equality 
of (4.9) and (4.10) follows the equivalence of a rotating pressure dipole 
and a vortex system. This system consists of a rotating bound vortex 
element of length dg stretched along the radius 9 and of trailing vortices 


Fig. 2. The rotating vortex element. 


along helicoidal lines (fig. 2). The first term in (4.10) originates from the 
bound vortex element which, in order to induce the correct velocities, 
must have the strength 

5 
Sore Y= RO Pay de 
The second term is caused by the two trailing vortices with strength + y, 
which are close to each other. Their resulting induced velocity can be 
obtained by differentiating the induced velocity of the outmost trailing 
vortex with respect to o. 

When the velocity fields of both representations (pressure dipole (4.9) 
and bound vortex with trailing vortices (4.10) are equal, then also the 
external forces exerted by these entities must be the same. We know 
that the external force exerted by the pressure dipole (4.3) has the magni- 
tude 8. From this and from (4.11) it follows that the force 7' on the rotating 


bound vortex with intensity y is 
(4.12) T = yuU (1+ ao?) do. 


The quantity U(1+a%0?)* can be interpreted as the local velocity of 


the bound vortex. 
Equation (4.12) states the law of Kutta-Joukowski for a rotating 


bound vortex element. 


5. The integral equation for the distribution of pressure dipoles 
It is our aim to formulate some integral equations which have to be 
satisfied at the blades of the screw. First we will derive such an equation 
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from the point of view of the pressure dipole layer. Then we have to 
introduce the velocities induced by a pressure dipole (4.9) into condition 
(2.6). The strength S of the dipole becomes a function of the coordinates 
£ and o. After integrating over the blade we obtain 
of 
roel 
4c epU?r = = 
3 *—~* -3¢ aor —r sin (a(t —x) + —)} {art —o sin (a(t—x) +¢)} 
= lim {1 5, 0) f [ape gmp comtale—aleone 


par Oo 


te ; {a?ro + cos (a(t — 2x) + ¢)} 1, def 
a ~ aE gi oon tea) P| See = 


a—é : 
lim ff S(E, 0) J F(x, r, yp. 0, t) dtd&do ae 
gar 0 —oo lim ff S(é, 0) K* (a, r,Q; E, 0) dEdo. 


paz O 


The limit mg — ax means that the point a, r, m, in which we consider the 
induced velocities, tends to the helicoidal surface H — 0. When we simply 
put the value g equal to ax into K*(x,r, y, §, 0) we find 


( K(2, 1, £, 0) = K*(a, 1,02, €, 0) = 


e—€ 
(5.2) r. [ ere (art — 0 sin at) (a2ro +- cos at) ] 


~ {724724 @2—2rocosar}s {72 + 72+ 9? — 2ro cosar}* 
—0o 

It can be seen that for §<«# and 9 — r this function becomes infinite. The 
singularity arises from a small part of the range of integration in the 
neighbourhood of t=0. For the study of this singular behaviour the 
range of integration in (5.2) may be changed into —«#<t<a, where « 
is some sufficiently small but fixed positive number. Putting 9=r+y 
where » is a small quantity it is valid to expand the numerators with 
respect to t and ». We obtain integrals of the form 


+o 
: 1 rmdr — — a 
=) 4 | {v2 72+ (r+)? — 2r(r+) cosar}a2 ~ <= 2 oUF) 
+o 
- ae TAT - 
(54) = | {7?(1+ ka®r2) + y2}0/2° 0<k<1, 


— 


where / and m are integers. It can easily be shown that for l+m>q-—1 
this expression remains finite when » — 0. Hence in (5.2) we have to 
expand the numerators of the first and the second fraction up to and 
including terms of the third and first order respectively. By doing this 
we find, after an estimation of the resulting integrals the following 
first two terms of K(a, r, é, 0) 


- S : i Sib aes ra* 
(5.5) K(a, r, é, 0) roo (eee) F<a,@>r. 
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This singularity however cannot be integrated and hence we have to 
follow a different method for giving a meaning to the limit stated in 
(5.1). This will be discussed in the following. 

We divide the area of the blade into three regions (fig. 3). The strip 
QO, defined by xi(r)<é<ax—y, lo—r|<f, the rectangle O2 defined by 
|v —&|<y, |o—r|<B and the remaining part 03. In this definition y and f 
are arbitrarily small positive quantities. 


MUMME EL a WW. 
CN AN 


r(r) © é 


Fig. 3. The domains of integration. 


The domain O03 does not yield any difficulty in the limiting process 
gy — ax because the point £, g remains at a finite distance of the point a, r. 

The integration over 0; when first » aw and then £0 can be 
written as 


(5.6) lim ti SE 0) L Fle. r. 9, 0.0) dr} dd 


P—-aL x (r) r—B 


where « is again a sufficiently small but fixed quantity. This is valid since 
only the singular part of the kernel can yield a contribution to the integral 
over the strip of vanishing width, 6 — 0. Then it is possible to expand 
the function S(&, 0) and the numerators of F(x, 7, p, 0,7) with respect 
to the small quantities 


(5.7) 6=y—at, Y¥=o-T, T. 


The integrals which appear are of the type (5.3), only now a factor 6” 
arises in the numerator. They can easily be estimated. It turns out that 
again a finite number of terms need to be considered which for the 
limiting case, first 6 > 0 and then f — 0, yield a contribution. After this 
it is possible to expand the denominators where also a finite number of 
terms need to be taken into account. In this way we find after a complicated 
but elementary calculation for the contribution of (5.6) 


(5.8) sper | Dien) de, p> 0. 


a, (T) 
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The integration over the small rectangle O2 can be treated in a similar 
way. Its contribution can be written as 
(5.9) lim aie Tse, 0) { "| Fle, r, p, 0, t) dt} d&do. 

pan x—y r—B = 

The lower boundary —ce for the integration over 7 is replaced by the 
fixed value —y<0 because the range of integration —co<t<—y does 
not yield a finite contribution when f —> 0. We now expand S(&, @) by 
a Taylor expansion for two variables and integrate partially with respect 
to €. Then we arrive at similar formulae as we found for the region QO}. 
It turns out that the contribution of the integral over Oz becomes 


(5.10) a [Sends Bo. 


z-y 


This is the extension of the range of integration in (5.8) from x—y to x. 
Using (5.8) and (5.10) we obtain for (5.1) 


r, ~'o) 
[ovr =m het\t 


r; r+B a0) 


-S(E, 0) K(w.1, & 0) dédg +E EO" | sie, nde], 


2, (r) 


(5.11) 


where K(x, 7, €, @) is defined in (5.2) and its singular behaviour is given 
in (5.5). This is called the Hadamard principal value of the integral. 

For the case a=0 we have no rotation of the blade and hence our 
formulae must reduce to those of the lifting surface theory for wings. 
The coordinates x and r, resp. € and o on the blade of the screw then 
become cartesian coordinates. From (5.2) it follows that (5.11) reduces to 


r—-B omy CO) 
[suo =n [J a pil 


r; r+B a0) 


(5.12) 


—~—, 
— 


hes \ déd 4 sen ) de]. 


{(w— &)? + (r—)?}*J 


_ SUE @) 
(r—@)* 


x (r) 


where the function ®=r(d//dx) is the angle of attack of a surface element 
of the wing. This formula agrees with the basic equation in [1]. 


6. The integral equation for the distribution of bound vortices. 
Here we start from the second formulation of the induced velocities 
(4.10). These have to be substituted into equation (2.6). After performing 
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a partial integration with respect to 9 we obtain 


Axe U2r ie = lim S(E, ae er Mi 8) a el sede) 
| ow p> ax [i ‘s @) {(x €)2 + 72 aio “— 2re cos (p— ag) )}*ls ac40: 


R oS(é, -s . é 
(6.1) { — pp °S&2) are, », 9,8, 9) dédo+ f (SE, e(6)) Ma, 1, o &, eu(§)) + 


) g 5 


v 


— S(E, ea(€)) Mle, r, @, & on(6))} ds] 


( M(x, 7, p, €,@) = 

\ (a—€&) 

| Be | {(i — a?r?)o cos (a(t — x) + ¢) +-agr sin (a(t — 2) +g) —r(1 — a?o?)} d 
{t?-+ 1? + 9? — 2ro cos (a(t — a) + p)}*/2 e 


Yq and x, denote the largest and smallest value of the € coordinate and 
ou(€) and oa(é) describe the circumference as function of & (fig. 4). 


Van LY, Ws WN 
+f CL, bam 


Py — 


When we consider the case p= az the kernels in (6.1) possess singularities 
which can easily be estimated. We find for the kernel of the first integral 


ar? sin (p —a€é) + (x— &) cos (p —aé) ae 
\ im E {(a— 6)? + 72+ 9? — 2ro cos (g— Zi lowe se 


(6.3) or . 
ee ee) ee 
TE or )\(E—a)* 4 (9 — 1) ye 


The singular behaviour of the function M(x, 7, y, €, 0) is given by 


+ 2(1+ a2r?)"/s 
\e=7) 


(6.4) lim W(x, r, ax,&, 0) ; Exe 0, 


oT 

By excluding the strip |r—o|</ from the domain of integration, it is 
possible to carry out the integrations in (6.1). 

From the formulae (6.3) and (6.4) for the singularities it follows that the 

limit ~ + 0 exists. However when the integrals can be interpreted in 

this way, we are not sure that the result is correct. There remains a 


21 Series B 
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possibility of contributions of functions which for p=axv are zero every- 
where except for g=r (6 functions of Dirac or its derivatives). That this 
does not happen follows from a more carefull passing to the limit, first 
g > ax and then p — 0. . 

The last integral in (6.1) is one dimensional with ¢ as the variable of 
integration. It is more natural to introduce ¢ instead of € because then 
for ri<o <70 (fig. 4) we cover the whole leading edge and the whole trailing 
edge separately. 

Ultimately we come to the following form for equation (6.1). 

p< §, (0) 
J 
£,(0) 


4reuU?r S = lim (f+ 


To 
. 


po ‘ a 
far? si == £) xf = 
(5, @) {(a — &)2 + 72+ 92? — 2rp cos a(x — €)}** dédg + 


r—B oy Fle) 


—({ + [) i} dS(E2) (x, r, aa, E, 0) dédo + 


(6.5) oo 
T; r+B &(o) 
Nihay es To ap ) 
1 el + [) S(E,(o), 0) M(x, r, ax, &(0). @) ( ae ) do + 
T; r+p 


To 


= | S(E(@), 0) M(x, r, ax, &,(0), @) (=e) do | 


% 


We shall now discuss shortly the well known physical meaning of the 
different parts of the right hand side of equation (6.5). The first integral 
represents the velocities induced by the bound vortices of which in fig. 4 
a typical one is drawn, the thick line AB. The second integral provides 
the velocities induced by free vortices, lying on screw lines, caused by 
variations in the o direction of the bound vortices. In the figure these 
free vortices are drawn as horizontal lines starting at AB and stretching 
to the right. The last two integrations in (6.5) give the velocities induced 
by the free vortices which arise by the ending of the bound vortices at 
the circumference, point A or B. The strength of such a free vortex, 
denoted in fig. 4 by a thick horizontal line, is equal to the strength of the 
bound vortex at the circumference. The relative rotation with respect 
to the bound vortex is indicated by arrows. 


7. Final remarks 


In lifting surface theory there are two main types of problems. First, 
the problem of determining the shape of the surface for a given load 
distribution. Second, the inverse problem where the surface is given and 
it is asked to determine the pressures exerted by the fluid. In contrast 
to the lifting surface theory for airplane wings, the first problem is more 
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interesting from the point of view of ship screws. Mathematically it 
means that in our integral equations (5.11) and (6.5) the function S(E, 0) 
is known. In order to calculate the function /(x,r) which determines the 
angle of attack and the camber of the blade sections, we have to carry 
out the indicated integrations. This is rather cumbersome owing to the 
complicated nature of the kernels. Equation (5.11) looks more simple 
than (6.5). However in this last equation the kernel functions involved 
are integrations over velocities induced by bound vortices along straight 
lines and free vortices on screw lines. These induced velocities have been 
studied rather intensively (induction factors) [5]. Hence it seems recom- 
mendable for numerical calculations to start from equation (6.5). 

Another point which has to be discussed is the extension of the theory 
to more bladed screws. Here we have to consider also the velocities 
induced by the other blades. We assume that the pressure distributions 
on the N equally spaced blades are identical. Then we can simply replace 
(5.1) for instance by 


d . Bete. 2: 
(7.1) 4neU2r a — lim ff S(é, 0) y K(e rp+—.é, ) dédo. 


pax 0 


The added part of the kernel has no singularities, because the points of 
one blade are at a finite distance of the points of the other blades. From 
this it follows that all our limit considerations and statements about 
integrability remain valid. The same can be done with equation (6.1). 

Summarizing the results of this paper, we can state that we have 
derived, for the one blade propeller, two equivalent expressions (5.11) 
and (6.5) which give the relation ship between the pressure distribution 
and the geometrical form of the propeller blade. It is easy to modify 
these equations so that they become valid for the screw with more blades. 
In general it may be expected that for an accurate numerical treatment 
of the two problems, stated in the beginning of this paragraph, electronic 
computers will be essential. 


Netherlands Ship Model Basin, 
Wageningen 
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PHYSICAL CHEMISTRY 


A COMPARISON OF AN EMPIRICAL AND OF A KINETIC 
EQUATION, BOTH DESCRIBING A SELECTIVE 
CATALYTIC PROCESS 


BY 


R. J. A. M. VAN DER BORG 


(Communicated by Prof. J. H. pe Borr at the meeting of September 26, 1959) 


1. INTRODUCTION 


Recently a method for simplifying process development has been 
proposed by WEBER [1]. Analysing a great number of selective chemical, 
mostly catalytic, processes he discovered a graphical statistical method 
by means of which all the examined processes could be described. 

In this method an isosceles rectangular diagram is used, in which the 
result of each individual experiment is indicated as an experimental point, 
(x, y). in such a way that the fraction (y) of the starting material A which 
is transformed into the desired product B is plotted vertically and the 
fraction (x) of A which is converted into an undesired by-product C is 
plotted horizontally. WEBER discovered that all points (x, y) of a series 
of experiments, in which only one reaction parameter was varied, as e.g. 
the contact time, could be described by an equation: 


a(1—z2) 


(1 Ears 


a and 6 being empirical constants. As, apparently, all series of practical 
experiments may be described by this equation, each series giving one 
set of the empirical constants a and 6, the number of experiments of each 
series may be drastically reduced, as is most clearly brought forward by 
WATERMAN [2]. This graphical tracing method, together with the empirical 
equation (1), therefore, makes it possible to determine the optimum 
process conditions after far less experimental observations than before. 


2. COMPARISON OF THE EMPIRICAL EQUATION WITH A KINETIC APPROACH 


In discussing various applications WarTrERMAN [2] and WeBeER [1] 
compare, in some cases, equation (1) with equations which are derived 
from kinetic considerations. The main example they give is the catalytic 
hydrogenation of linoleic esters, which reaction has been investigated, in 
earlier years [3], in Waterman’s laboratory, and which can be described 
as a consecutive first order reaction [4]: 


k *, stearic ester (S). 


linoleic ester (Z) “> oleic ester (OQ) - 
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In a selective hydrogenation oleic ester is the desired product B and 
stearic ester the undesired by-product C. In discussing various types of 
selectivity WuHrELER [5], indicating this type as type III selectivity, 
derives, from kinetic considerations the following equation: 


S 
(2) wp = yoy (1— aa) [(1— oa) 0/9) — 1] 


where xg is the fraction of initial A converted into B, «4 is the fraction 
of A reacted and S is the so-called selectivity factor k1/ke. 

Applying both equations (1) and (2) to their hydrogenation experiments 
of linoleic ester, WEBER [1] and WaTreRMAN [2] derive the conclusion 
that the hydrogenation experiments at high temperatures (180° C) can 
be described by both equations. The hydrogenation experiments at low 
temperature (20° C-90° C), may also well be covered by the empirical 
equation (1), equation (2), however, fails in this case. 


3. A NEW APPROACH 

Equation (2) as derived by Wheeler contains only one reaction constant, 
S. Equation (1) is characterised by two empirical constants, a and b, 
A slight extension of the kinetic scheme, underlying (2), may introduce 
a second kinetic reaction constant. 

Instead of this scheme: 


(A) L“+0-"S8 


another scheme will be used, allowing a direct conversion of L into S: 


B Ie 
a Lo——>8 


If all reactions may be described as first order reactions, the following 
rate equations hold: 


(3) dL [dt = —ki L—kgL 

(4) dO/dt = ki L—k,O 

(5) dS/dt = kgO + k3L 
whereas, when weight or mole fractions are used: 
(6) L+0+4+8 =1 

(7) O=-y 

(8) S= 2. 

Initially 

(9) Leal 


and _ finally 
(10) S=1, 
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As shown in the mathematical derivations at the end of this paper, 
this set of equations leads to: 


(11) (l—x—y)! = 1-x—my 
where 
12 Ses 
( ) l 1+ks3/ki 
and 
(13) m = kelky —ks/ky. 


ey 


4, EVALUATION OF EQUATION (11) 

The constants J and m may be evaluated in the following way. The 
experimental data are plotted in a rectangular diagram, as indicated in 
the introduction and, as a first approximation, a curve is drawn in the 
way it is expected to run. The values of x and y in the maximum (see 
fig. 1) of the (real) curve give k2/k; (mathematical derivations): 


l= = Yrna: 
ay eer re 
The direction of the curve in the origin (mathematical derivations) gives 
kg/k1: 

(15) kg/ky = cotg «. 


The provisional curve gives provisional values for | and m. The exact 
values may be estimated from two experimental points (7, yi) and 
(a2, y2) as follows. Elimination of m in equation (11) gives: 


(16) ay a2 — yz)! me a) +(1—a)). 


Both the left hand side of (16) and the right hand side of this equation 
are then drawn as functions of / (values of / round the approximate one). 
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The intersection of both curves (see fig. 2), f(J)iere and f(U)rignt respectively, 
gives the exact value of /, whereupon that of m can be calculated. 


eth) 


5. SOME EXAMPLES 


Some of the curves of low-temperature hydrogenation, which cannot 
be described by Wheeler’s equation (2), but which can be described by 
the empirical equation (1) are taken from the figures of WeBER [1] (his 
fig. 28) and WATERMAN [2] (his fig. 93) and shown in fig. 3. The constants 
a and b of the drawn lines are given in table I; they are taken from table 


VIII of Weezer [1] or from table XVII of WaTERMAN [2]. 


1.00 


o 2.5% Ni-guhr 90 °c, 1 atm 
o 10 &% Pt-carbon 30 °c, 1 atm 
e 10 &% Pt-carbon 20 °C, 100 atm 
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Table I also gives the values of / and m calculated along the lines of 


section 4 of the present article. 


TABLE I 
rr ETT 


| 


Seanad e | Pressure ‘Temp. | | 
ymbo Catalyst | aan eb OG) | b l m ke/ki | ka/ki 
Gr 2696. |) eee 
| Ni-guhr 1 90 | 0.117 | 0.712] 0.150) +0.049 | 0.17 | 0.12 
Er | 10.9% | | 
Pt-carbon 1 30 | 0.340 | 0.398 | 0.178 | —0.142 | 0.25 | 0.39 
| | | 
@ | 10 % | | 
| Pt-carbon 100 20 | 0.747 | 0.128 | 0.373 | —0.110 | 0.65 | 0.76 


With the aid of the values / and m fig. 4 is drawn. Moreover ke/k, and 
k3/ky can be calculated; these values are also given in table I. 


6. COMPARISON BETWEEN EQUATIONS (1) AND (11) 

Figs. 2 and 4 show great similarity. An advantage of the curves derived 
from (11) is that they remain inside the triangle. (1) describes the obser- 
vations only in so far as the curve remains inside the triangle, the hypo- 
thenuse being the continuation of the curve. 

A comparison of the numerical values derived from (11) and (1) is 
shown in table II. 

TABLE II 


2 Oo ®@ 
C y pO] y y y y 
(11) dq) | dy (1) (11) (1) 
| | 

0.062 0.362 | 0.362 — =e ee. 
0.100 | 0.486 | 0.479 | 0.226 | 0.236 Oolde 0.118 

0.174 | — = im 0-352, | =0.352 | a = 
0.183 a =e 0.194 0.194 
0.200 0.627 0.618 | 0.385 0.382 0.207 0.207 

0.277 0.637 0.637 ae ae = ae 
0.300 0.632 0.634 | 0.475 0.458 0.270 0.267 
0.400 0.576 0.597 | 0.495 0.481 0.306 0.301 

0.419 = = 0.493 0,480 a = 
0.500 0.493 (0.500) 0.460 0.464 0.314 0.308 
0.600 0.398 (0.400) 0.388 | (0.400) | 0.294 0,292 
0.671 ——- — _ = 0.265 0.265 
0.700 0.300 | (0.300) 0.298 (0.300) 0.250 0.250 
0.850 0.150 (0.150) | 0.150 | (0.150) | 0.142 0.149 


Where (1) falls outside the figure the values of the hypothenuse are 


given in brackets. . 
The y-values have been calculated from the x-values with the help of the 


parameters given in table I. 
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Despite the similarity between the two series of numerical values we 
have not succeeded in reducing (1) and (11) to each other. 


7. MATHEMATICAL DERIVATIONS 
a) Derivation of equation (11) 


From the equations 


(7) G=¥ 
(8) S = 2 
it follows: 
L=1-2-y. 
We put 
(8a) l-w=u. 


For equation (3) we can write: 


d : 
(17) WY) = _ky(u—y) —ka(u—y) 


equation (4) becomes: 
(18) dy/dt = ky(u—y) —key 
and as da/dt= —du/dt, we obtain from (5): 


(19) du/dt = —koy—ks(u—y). 


Division of (18) by (19) gives: 


(20) dy = kyu —ky = koy 
du —ksu — key +ksy 

Suppose 

(21) y = vu 

hence 


dy/du = u(dv/du) + v. 
Substituting this into (20) we obtain: 


dv = ky —kiv —kov +khgv + kev? — kav 
du —k3 —kov +k3v 


The variables can now be separated and the equation be integrated : 


(22 —ks+-(ks —ka)v ach du 
} — (ks —ke)v® —(ki +k —ks)v +k dv = u- 
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The left hand side of (22) can be simplified and be split up into partial 
fractions, leading to: 


ke [ dv _ ky tks f dv < r du 
ki—ko+k3 J v—1 ee v+ki/(k3—ke) | Uw 


ke kitks 
ki —ke +k lnjjo= i es rere |v+ ki/(ks —ke)| = nu + In Py 


Eilat Haan Py y) h1-kotkg) | (ky+k 3) (1—ko/ (ky +h w 
— vi1—Kot+k;: Y = DY —K2 tk: = ——— 
|v + k1/(k3 —ke) | raya sat at Pat ae Pe wkal (ky +s) 


hence 


ac —w| Wa! (1 +s) 


23 Fila al al 
(23) u|v + ki/(k3 —ke)| 


wi 


Substituting (8a) and (21) into (23) we obtain: 


ly—1 +a] Fa! (es +s) me P 
ly-+(1 —a) ki/(ks —ke)| Ze 


Because the curve must remain inside the triangle, the equation can be 
written: 


(1 —a—y)*2/(eitks) = Psll —x—y(ko—kg)/ky]. 


The initial condition (x; y)=(0; 0) gives P3=1. 
Substituting (12) and (13) we finally obtain: 


(11) (l—x—y)! = 1—x—my. 
b) Another derivation of equation (11) by H. A. G. CHERMIN 
Equation (17) immediately leads to: 
(24) u—y = exp {—(ki+hs)t}. 
Substituting (24) into (18) we obtain: 
dy/dt+key =k, exp {—(ki+ks)t} 
which leads to 


t t t 
y= exp (—f kedt)[C+ J kiexp {—(k1+ ks) o} exp ({ kodt) dt]. 
0 0 0 
When t=o, y=o, hence C=o, hence 
t 
y = exp (—het) f kiexp {(k2—hki — ks) t} dt 
0 


y = exp (— kat) aS [exp {(k2— hi — ks) t}— 1] 


(25) las ao [exp {— (ki + ks)t} — exp (—hel)]. 
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Putting 
(26) (ky +ks)t = 
and using equations (12) and (13) for / and m, we obtain 
(27) y = 1/(m—1) {exp (—t) —exp (—/t). 
(24) together with (26) gives: 
(28) 1—x—y = exp (—7). 
Substituting this into (27) we obtain 


(m—1)y = (1—a—y)-—exp (—It) 


ee exp (—Ilt)=1—2—my. 


As 
exp (—I/1) = {exp (—7)}! 
we obtain from (28) and (29): 
(11) (l—x2—y)! = 1—x—my. 
This derivation has the advantage that also y and (w—y) are obtained 


as functions of ¢ (equations (24) and (25) respectively). 


c) Derivation of equation (14) 
In the maximum of the curve the first derivative is equal to zero. 
Division of (4) by (5) gives 


dO kxL—kO 


(30) dS ksO+ksL 


hence 
Ky Lynax — K2Omax = 0 


ka 1 —2max — Ymax 
14 ka 1~®max ~Ymax_ 
( ) ky Ymax 


d) Derivation of equation (15) 
Equation (15) likewise follows from (30). For L=1 and O=0, we obtain 


dO/dS = ky/kg = tg « 
(15) k3/k, = cotg «. 


The author wishes to thank Prof. J. H. pz Borer and Prof. G. MEYER 
for their valuable stimulation and comments, and especially also Mr. H. 
A. G. CuEermin for his contribution towards the mathematical derivations. 
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PHYSICAL CHEMISTRY 


THE KINETICS OF THE SELECTIVE FORMATION OF THE 
INTERMEDIATE PRODUCT IN TWO CONSECUTIVE REACTIONS 


BY 


J. H. DE BOER anp R. J. A. M. VAN DER BORG 


(Communicated at the meeting of September 26, 1959) 


1, INTRODUCTION 
In a previous article [1], one of us (vAN DER Bora) has shown that 
the formulation of a selective catalytic process: 


(x) As FR 


where B is a desired intermediate, which can react further to give an 
undesired by-product C, may substantially be improved when also a 
direct transformation of A into ( is taken into consideration: 


a B., 


4 k 3 


It could be proved that a kinetic equation based on scheme (f) describes 
the experimental data equally well as an empirical equation which is 
recommended for its universal use [2]. In such a case the equation which 
is based on kinetic considerations is to be recommended in favour of the 
empirical one, because its constants have a physical significance. 

In the present article the physical meaning of schemes («) and (#) and 
of the equations, which are based on these schemes will be discussed. 


2. EQuaTIONS 

a. The general equation 

If the three all-over reactions, for practical purposes, can be described 
as first order reactions, scheme f leads [1] to the following equation: 
(1) (1-a—y)! = 1-a—my 


where 2 is the fraction of A which is converted into C and y is the fraction 
which is converted into B and where: 


ke/ky 

2 ee I 
(2) Be ee = 
and 

(3) m= ko/ky = k3/ky. 


The parameters of equation (1), therefore, are the relations between the 
rate constants, k2/k; and k3/k; respectively. 
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b. The special case of k3=0 


If k3 + 0, scheme (f) reduces to scheme (x), which is, therefore, a 
special case of the more general scheme ({). Substituting /3—=0 in equations 
(2) and (3), we obtain: 

b= m = kolki, 


and equation (1) reduces to: 
(4) (1—a—y)W = loa By 
Bi 


Equation (4) is identical with the equation which WHEELER [3] derived 
from the kinetics of scheme (a), viz.: 


(5) kp = os (=e) [(d—es)0-9/9) 1) 


where ag is the fraction of A, converted into B, 
«4 18 the fraction of A which has reacted, 
and S is the so-called selectivity factor; S =ky/ke. 


The identity of (4) and (5) can easily be seen, when we note that: 


op=y and o4g=a+y. 


c. The special case of k2=0 
If k2=0, we obtain: 


f=, 
os kg/k1, 
and equation (1) reduces to: 
(6) ccameag onc 


Scheme (f), in such a case, reduces to scheme (y) 


(y) A 


Wheeler’s equation (7), derived [3] from scheme (y), viz.: 
(7) op = Sxc, 


where: S=ky/ks, ; 
is identical with (6), when we note that in this case 


op=y and wc= a. 


3. THE PHYSICAL MEANING OF SCHEME (/) 
Suppose A to be a molecule, which, on a heterogeneous catalyst, can 
react with a molecule M to form a product B, whilst this reaction may 
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be followed by another reaction with M, on the same catalyst, to give 
product C'. To give an example, A may be an organic aromatic compound, 
such as benzene or naphthalene, and M may be oxygen. On a vanadium 
oxide catalyst product B may then be maleic anhydride or phthalic 
anhydride, respectively and C is carbon dioxyde. In a hydrogenation 
process, A may be an acetylenic compound, B the corresponding ethylenic 
compound and C the saturated compound, or, if A is linoleic ester, B 
may be oleic ester and C stearic ester. In all these hydrogenation processes, 
M is hydrogen and the catalyst is e.g. nickel or platinum in a suitable form, 
In order to execute the partial process A — B molecules of A and of 
M have to move to the surface of the catalyst (diffusion to the outer 
surface and, if necessary, diffusion into the pores), they have to be 
adsorbed on that surface, whereupon they can react to form adsorbed B, 
which has to be desorbed, followed by a diffusion of B from the surface 
of the catalyst to the bulk of the homogeneous (liquid or gaseous) phase. 
Apart from all diffusion processes and apart from the chemisorption 
process of M, the partial process A — B involves: 
a) chemisorption of A, 
b) reaction of chemisorbed A and M to chemisorbed B, 
c) desorption of B. 


The partial process B-—->C, treated similarly, involves: 
’ J? 


d) chemisorption of B, 
e) reaction of chemisorbed B and M to chemisorbed (C, 
f) desorption of C. 


The direct route A —>C of scheme (8), however, comprises: 


a) chemisorption of A, 

b) reaction of chemisorbed A and M to chemisorbed B, 
e) reaction of chemisorbed B and M to chemisorbed ©, 
f) desorption of C. 


We assume that the chemisorbed state of B, produced in reaction b) 
is identical with that which results from the chemisorption process d); 
in that case step e) in the partial process B —>C is the same as in the 
direct route A > C. 

The essential difference between A —- B->C and A -—->(C, therefore. 
is the behaviour of the intermediate B. 

In the latter, direct, route, the intermediate B is not desorbed, and, 
consequently, does not contribute to the content of B in the bulk phase. 
If, during its time of adsorption [4a] a molecule M can react with the 
adsorbed intermediate, step e) is performed. If, however, its time of 
adsorption is so short that step e) does not take place, B is desorbed. 
The time of adsorption depends on the activation energy for desorption 
and on the degree of mobility of the adsorbed species [4b]. The rate of 
reaction of the intermediate depends on similar entities. 
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4. ‘THE FATE OF THE ADSORBED MOLECULE B 
In order to get an idea about the fate of the chemisorbed molecule B 


after step b) is performed, we shall apply some equations of the theory 
of absolute reaction rates [5] to this problem. 


a. Localized adsorption of B and M 

If we suppose both B and M to be adsorbed on definite adsorption sites, 
without translatory movement across the surface and without rotational 
freedom, the rate, v-, of step c) may be given as follows: 


OE Fs. 
(8) Po AO Ge 
B 


e-Lakr 
h 


where Cz is the surface concentration of B, F., is the partition function 
of the activated complex for desorption per cm?, Fs is the partition 
function of the adsorbed molecules B per cm2, Hq is the activation energy 
for desorption (per molecule), & Boltzmann’s constant, h Planck’s constant 
and 7’ the absolute temperature. 
Similarly the rate, v., of step e) following b) in the direct route A + C 
is given by: 
on Cem, kT _ (Fs)am 


(9) Ve ‘7 "er x ae 6 Fact/kT 


where Cy, is the surface concentration of MW, L is the total number of 
adsorption sites per cm?, s the number of direct neighbour sites, sur- 
rounding any adsorption site, (F.)gy4 the partition function of the 
activated complex for reaction, per cm?, F'gy the partition function for 
the direct neighbouring pairs of B and M per cm? and Kact the activation 
energy for reaction. 

Dividing (8) and (9), we obtain: 


vw _ Bs y Psu _ —(Eg-Eact) /kT 
(10) fe Ogee ZW Prd ac 
Cw can be written as Ly, where 6 is the degree of occupation of M. 
Because of the localised character of the adsorption, #g=1 and Fgy=1. 
F.. and (F.)gm need not be unity, but we assume them to be roughly 


equal, in which case (10) simplifies to 


(1 1) dt = —_ Cm (Eqg-Lact) /kT 


Ve Ms 


b. Mobile adsorption of B and M 

If we suppose B and M to be in a mobile adsorbed state, we obtain for 
the rate ve the same equation (8) as in the immobile case. However, 
F. and Fz are not equal to unity in this case. 

For v,, we obtain: 


Nig ean 0 aes ae ree 
(12) Ve = CpXCuX 7 x Pasian << e- Eact [kT 


22 Series B 
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where (F.,), is the partition function of the mobile activated complex per 
cm2 and Fz and Fy the partition function of B and M, respectively, 
both per cm?. 

Dividing ve and ve, we obtain: 


(13) Ve 1 By X Pu ong 


aoe, peat Mee Eact) (kT 
Ve Om (fs )r 


For Cy we write again L6,, and we assume F=(F',)e- 
Equation (13) then simplifies to: 


(14) = = ne e~ (Ea-Eact) /kT 

c. Comparison between equations (11) and (14) 

The pre-exponential factor in equation (11) contains, apart from yy, 
the number s, which may be 3, 4 or 6. The pre-exponential factor in 
equation (14), however, contains, apart from §, the proportion F'y/L. 
L is of the order of magnitude of 101° em~2, whilst: 


2amukT 
(15) Fu= a x (Fa) rot X (F'm) vin 


where my is the mass of a molecule M and (Fyy)rot and (F'y4)yin are the 
partition function of the rotational and vibrational degrees of freedom, 
respectively. 

If the molecular weight of M=20 and 7’'=300° K, (15) gives: 


Fy = 20x 1016 (F'yy)rot (F'at)viv 


and the pre-exponential factor of (14) may be roughly 108 times the pre- 
exponential factor of (11). Even if 1/ =H (hydrogen atom), as may well 
be the case in hydrogenation reactions, Fy—10' at 300° K and 
Fy=1.5x 106 at 450° K (rotational and vibrational contributions can 
be neglected) and the pre-exponential factor is still 40 to 60 times higher 
in the case of mobile adsorption. 


5, INFLUENCE ON SELECTIVITY 


If the intermediate product B is the desired product and C is an 
undesired by-product, k; should be large and, both, ke and kg should be 
low. A low value of kg is obtained if v, is high with respect to v. 

If, in equation (11) or (14), Ha=Hact, the proportion v/v, is only 
governed by the pre-exponential factor of both equations. In such a case 
localized (immobile) adsorption would hardly give any selectivity, but 
mobile adsorption would result in an appreciably high degree of selectivity. 
Mobile adsorption is promoted by raising the temperature and by 
choosing a catalyst not adsorbing too strongly. 

It is a remarkable result of these considerations that raising the 
temperature should promote selectivity. This, of course, would be more 
pronounced if Ha < Fact. If, therefore, the activation energy for desorption 
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of the intermediate B is smaller than the activation energy of the further 
reaction of this intermediate, selectivity is certainly promoted by raising 
the temperature. At the same time, the condition Ha< Eact points to the 
selective character of a less active catalyst. 

Both equations (11) and (14) contain 6y7 in the denominator. Any factor 
lowering 6j,, therefore, also promotes the selectivity. 
Selectivity by suppressing the rate constant £3 can, therefore, be obtained 
y: 
a) higher temperature, 
b) a less active catalyst, 
c) a lower concentration of the substance M on the catalyst. 


oF 


6. THE SELECTIVE HYDROGENATION OF LINOLEIC ESTER AS AN EXAMPLE 


Literature shows that the selectivity of this hydrogenation reaction to 
form oleic ester, and preferably not stearic acid, is promoted by the 
following factors: 

a) higher temperature. WATERMAN [2] showed that the high temperature 
hydrogenation (180° C) could be represented by equation (5), hence 
by scheme (x). It means that, at that temperature, reaction kg is 
practically suppressed ; 

b) a less active catalyst; Nz is a more selective catalyst than Pt; 

c) a lower concentration of hydrogen on the catalyst, which may e.g. 
be obtained by a lower hydrogen pressure. 

If we take the selectivity factors S=k,/k2 for higher temperatures 
from the figures of WATERMAN [2] and, if we denote S=ky/(k2+ 3) for 
lower temperatures, where kz is still significantly high and we calculate 
S from the figures, given in the previous publication [1], we obtain the 
following table, which illustrates the relations very well. 


TABLE I 
| ie Hydrogen. 
Catalyst ve | pressure Ss 
| ; atm. 
10% Pt-carbon | 30 1 1.55 
20 100 | 0.7 
180 1 6 
2% Ni-guhr 90 1 3.4 
ee l ~25 
| 180 100 ~ 2 


7. "THE SPECIAL CASE OF SCHEME (y) 


If 1) scheme (y) is a special case of scheme (jf), step BC is suppressed, 
which means that adsorption of B from the bulk phase does not lead to 


1) A scheme (y) can also be visualized if the two reactions are quite independant 
from each other. We only consider here the special case that the step A > C takes 
its course via absorbed B. 
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reaction. A poisoning of the catalyst in such a way that B cannot be 
adsorbed properly may lead to such conditions. If the other conditions 
are such that also k3, hence step A->C is suppressed, ideal selectivity 
will be found. It is reported in literature that sulphur-poisoning of a 
nickel-catalyst may suppress the hydrogenation of oleic ester completely [6]. 

Another method of suppressing step BC may be applied by adding 
an indifferent substance with an adsorbability in between those of A and 
of B. Step B-—> C will then be suppressed, but not step A — C. Other 
measures, such as choosing a higher temperature, or desactivating the 
catalyst, or a lower concentration of hydrogen (in hydrogenation reactions) 
may be applied to that end. It is not improbable that some very special 
catalysts, used in the pharmaceutical industry operate along these lines [6]. 


8. OTHER FACTORS 


It must finally be pointed out that the kinetic expressions, leading 
to equation (1) are derived, assuming that the all-over rate reactions can 
be described as first order reactions. This may in practice often be done. 

Diffusion rates are not taken into consideration. In liquid media the 
bulk diffusion may influence the results of practical experiments. When 
microporous catalysts are used diffusion in the pore system may influence 
selectivity problems seriously. The influence of these factors on the 
selectivity problems discussed in the present paper, will be the subject 
of later studies. 
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PHYSICS 


ON THE POWER TRANSFER BETWEEN PARAMAGNETIC SPINS 
AND CRYSTAL LATTICE. Ia 


BY 


B. BOLGER 


(Communicated by Prof. C. J. Gorrmr at the meeting of June 27, 1959) 


INTRODUCTION 


The establishment of equilibrium between a system of paramagnetic moments and 
the crystal lattice can be characterized by a spin-lattice relaxation time 71. This 
has been the object of many investigations and recently the interest and that for 
related problems has gained new impetus due to the application of paramagnetic 
materials in low noise microwave amplifiers, so called “‘solid state masers’’. 

In the usual method to determine the spin-lattice relaxation time an alternating 
field is applied parallel to the static magnetic field and the complex susceptibility 
measured as a function of the frequency. With the aid of the thermodynamical model 
of Casimir and Du Pr [1], t1 can then be found. A survey of the experimental 
data up to 1947 has been given by GorTeER [2]. 

Apart from a few exceptions, the theory of VAN VLECK [3] describes the experi- 
mental results for the magnetically diluted salts fairly well. For the concentrated 
salts, on which most of the early measurements were carried out, the agreement at 
liquid helium temperatures is miserable. 

With the advent of paramagnetic resonance new techniques have become available 
to measure the spin-lattice relaxation behaviour. The system of magnetic moments 
can be heated, owing to an inadequate power transfer to the lattice, by absorption 
of microwave power. The study of the magnetic temperature as a function of the 
absorbed power in the stationary state, yields a value of the power transfer constant 
closely related to 71. This method is called the saturation method. Data obtained 
by EscHENFELDER and WEIDNER [4] seemed to indicate a large discrepancy between 
the results obtained by the two methods. 

It has been the purpose of the investigations described in this work to study the 
spin-lattice relaxation behaviour in a number of paramagnetic salts by the saturation 
method. A comparison is made between the results obtained by theory and by the 
relaxation and the saturation method. For the work carried out in connection with 
the construction of a maser reference is made to [5] and [6]. Only the results of 
that work pertaining to relaxation properties has been used. 
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ON THE GENERAL THEORY OF PARAMAGNETIC RESONANCE 
AND RELAXATION 


[1.1] The energy levels of paramagnetic ions 


Paramagnetic resonance deals with the transitions caused by a magnetic 
high frequency field between the lowest group of energy levels of a para- 
magnetic ion. The splittings between these levels will first be discussed. 

The case of free ions is relatively simple. The fact, however, that in 
crystals the energy levels are displaced by interactions with the sur- 
roundings makes paramagnetic resonance a fruitful source of problems 
in solid state physics. 

In this thesis we shall restrict ourselves to crystals in which the effect 
of the crystal field is small compared to that of the Coulomb interactions 
of the electrons with the nucleus and with each other, so that the self 
consistent field method can be used to describe the ionic wave functions. 
In ionic crystals this condition is fulfilled to-a high degree. 

We shall only consider the elements belonging to the iron transition group. 


[1.1.1] The Hamiltonian used to describe a paramagnetic ion in 
a crystal lattice can be divided into: 

a. the orbital energy #o due to the coulomb energy of the electrons. 
This energy is mainly due to the influence of the nearest atoms. The 
direct effect of ions at a larger distance is negligible and their influence is 
mainly felt via coupling with the nearest neighbours. 

b. The spin orbit coupling “os will be taken to be (ef. [1]) Hos=AL-S 
where L and S are the total angular momentum vectors of orbit and spin. 

A tensor type of interaction between L and § has been considered by 
Pryox [2]. 


c. The Zeeman energy 
(1.01) HAH z= B(H +h) (L+ 2S) 


where # is the Bohr magneton, H the applied static field and h a radio 
frequency field. 


d. The dipole-dipole interaction between two ions i and j can be 
written as 


(1.02) Hg — MAH _ 9 (Meta) (uy Tu) 
ry? ry? 


where rjj is the radius vector between ions i and j with magnetic moments 
wi and py. 


e. The exchange interaction 


(1.03) FE op = 2S 4384S; 


a7) 


is due to the overlap of the wave functions of ion i and 7; in this expression 
superexchange may be included. In more general considerations one has 
to add anisotropic exchange terms. 

The interactions of the electronic spin with the nuclear spin moment 
have been omitted as they are generally small and are not of interest in 
the present work. 

The terms d. and e. determine the coherence effects between the spins 
and will be discussed later. 

As the interactions #0, Hos and #7 are generally far more important 
than the coupling of the spins with each other, it is possible to consider 
the spins of different ions as being separate entities in a first approximation. 


[1.1.2] The perturbation calculation starts with the largest 
term # 0. The crystal field around an ion is expanded in Legendre poly- 
nomials and one considers how the orbital degeneracy is removed by 
going successively to lower orders of field symmetry. This has been done 
by the methods of group theory by Brerun, MULLIKEN and JAHN [3] 
[4] [5] and we will use their notation of designating the levels. 

For the case that the lowest eigenvalue of #o is non-degenerate, the 
perturbation calculation can be carried out by a method given by 
Pryce [6], by chosing a representation for L, diagonalizing #o and the 
so-called spin Hamiltonian contains the operator S and the matrix elements 
of L in the specified representation, corresponding to the lowest multiplet. 
It is now possible to define an effective spin operator S’ so that 2S’+1 
is the degeneracy of the ground state. The spin Hamiltonian can be 
expressed in terms of the operator S’. For the case of a degenerate lowest 
orbital level reference is made to the literature [7] [8]. 

For ions with L=0 one has to go to a high order perturbation, but for 
other ions the second order is usually sufficient. In the derivation one 
makes use of a general theorem that L has no diagonal elements for a 
singlet orbital state. 

The spin Hamiltonian, when spin-spin interactions (d and e) are 
neglected, has the following form: (cf. [7]) 


| H 5 = 2BH-S + (24)-1(BH + AS)(9— 2)(6H + AS) = 


(1.04) —9pH-S+S-D-S—f?H AH. 


If the ground state orbital wave function is denoted by <o| and the higher 
ones by <n| then 


= 2A oe 
(1.05) G-2)=S gH Lin <n|L;|0> 1,4 =H, 9.0L Z. 


The factor g is a tensor and so are D and A. 
In the third part of (1.04), the Zeeman energy is given by the first 
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term, the second term gives the zero field splitting caused by the crystal 
field and the third term represents the temperature independent para- 
magnetism. The total spin Hamiltonian must of course show at least 
the same symmetry as the crystal field. One advantage of using this 
method is that the resonance properties of an ion can be described by a 
small number of constants. For a more detailed discussion of the spin 
Hamiltonian and the values of the constants found experimentally in the 
various substances, reference is made to the review articles of BLHANEY 
and Stevens [7] and of Bowsers and OWEN [8]. 
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Fig. 1.1. Energy levels of KsCr(CN)s as a function of H. 


--- H//«w-axis. -—----- H // y-axis. H // z-axis. 

[1.1.3] As an illustration we give in fig. 1.1 the energy levels of 
Cr++* in CrK3(CN)g¢ diluted with the Co salt as a function of the applied 
magnetic field H along the different crystal axes, as calculated from the 
spin Hamiltonian. Cr+++ has a S=8/y and in this substance #s has the 
following form 


(1.06) H 5 =GgPH-S + D[S2—1/55(8 + 1)]+ E(S22 —Sy2) 


where g is almost isotropic and equal to 1.99. The values of the other con- 
stants are D=0.083 cm! and H=0.011 cm-l. The magnetic axes are 
slightly rotated with respect to the crystal axes. There are two Cr+++ ions 
ina unit cell with the ac-plane as a symmetry plane, 
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To find the paramagnetic spectrum at 4000 MHz (or 0.133 em-) for 
instance, one has to adjust the field so as to fit this energy between two 
energy levels. To show the anisotropy, a rotation diagram (fig. 1.2) in 
the ac-plane for the frequencies 1420, 3850 and 8500 MHz, has been 
drawn for some of the transitions. This diagram represents the loci of 
the field vector for which resonance takes place as a function of the 
angle « of H with the a-axis. 

The transitions are labeled by the corresponding energy levels, numbered 
from the lowest. 

c-axis 
2 


kOe 


3 


5 
, | 
O d=caxi's 
Fig. 1.2. Rotation diagram of Ks3Cr(CN)s, H in ac-plane. Corresponding energy 


levels next to loci, lowest is 1. 
1420 MHz. £—.--— 3850 MHz. —-—-— 9400 MHz. 


[1.1.4] The terms of the Hamiltonian of the paramagnetic ion 
depending on the position of the ion in the lattice such as in Wo, through 
a perturbation calculation, and in #q and #,, through the ion distances 
riz in (1.02) and (1.03), may be developed into a power series of the 
deviations u, and u,; of ion & and 7 from their equilibrium positions. 
Limiting ourselves to terms of second order, this gives 


(1.07) KH =H o —- > H 1;,(Ux) = > He 25(UKUj)- 
k Ij 


For the first term W%o, due to the ions in their equilibrium positions, we 
take the spin Hamiltonian (1.04) with the spin-spin interactions added, 
the second and third terms are respectively linear and quadratic terms 
in the displacement vectors uz. They give rise to the coupling of the 
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spin system with the thermal motions of the lattice. Their influence is 
several orders of magnitude smaller than 4%» and they can be regarded 
as small perturbations. To a first approximation the spins can be regarded 
as being decoupled from the lattice. 


[1.2] Macroscopic equations 


[1.2.1] A classical description of nuclear and paramagnetic 
resonance has been given by Briocu [9] for spins without zero field 
splitting. Consider an assembly of NV spins, all with the same gyromagnetic 
ratio y=gf/n. 

If the spin system is in thermal equilibrium with the lattice at a 
temperature 7’, the magnetization M will have the value given by the 
Curie law, ‘ > 

OH — Ny*h2S(S+1 
(1.08) M => = abcd 
where H is the applied magnetic field. 
The equation of motion of M without interactions with the surroundings is 


dM 
ae y(M x H) 

which expresses the fact that M precesses around H with the Larmor 
frequency : 

* y w= yH. 

Buiocw introduced damping terms in the equation of motion of M. 
Due to the interactions the different electrons will precess at slightly 
different frequencies. The components of M perpendicular to H will 
therefore get out of phase after a time ts, the so-called transverse 
relaxation time. Then 

d\|Mz y| |\Mz v| 
¢ al ere , 
8) dt T2 
where H is taken along the z-axis. 

The relaxation time 1, of the component of M parallel to H will in 
general have a larger value than t2, as a change in the z-component 
entails a change of the Zeeman energy and this can only happen through 
a coupling with the surroundings. 


(1.10) ¢M,_ _ Ms—Mo 


a TL 


We shall now write down the equations of motion when a linearly 
polarized radio frequency field h perpendicular to H is applied to the sample, 


dd Mz, M 
ae = y[M xX Hey = “ne 
(1.11) 2 
dM; 


= M;—Mo 
ar — VIMxH], — Mae 


with Hz=h/2 (e+e); Hy=0; H,=H. 
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The steady state solution has the following form 


4? +75? [ wp? + 1/t22 —w2]? 
— Miyh[io +1/r2] eit 
Y  (t@ 1/2)? +22 


M,= My? Ahre?| (wz? + 1/t22 —w?) — (2iw/t2) Jeiot 
x 


tars 4m? +72? (wp? + 1/T22 —w2)2 
ay 4° +72?( wz? + 1/t22 —ao?)? +1/ay2h2r172(wo? + 1/t22 +02) 


M, 


where the high frequency motions of M, have been eliminated by averaging 
over many Larmor periods. 

Near the resonance frequency w?=y2-+ 1/122, slightly displaced because 
of the damping, we have 

; M ht2(At2 +7) eiot 
(1.13) Niet aa ee 
2 1+ (Ate)? +1/ay2h2t172 
where 
woL+1/t22-—w2=2wA. 


Introducing the complex susceptibility y= y'’—iy'’ so that M= 7H, and 
relating this to the power P absorbed by the system from the r.f. field, 
one obtains 


(1.14) P=1/omy"h?2 =1/2mh-Im M. 


From (1.13) we find 


NN (fiw) 2y2h?r2.S' (S41) 

o>) on 12kT[ 1+ (At)? +1 /ay?2h2r 02] 

Upon increasing h we have initially P « h? as expected, but for values of 
1 


h?y——_—— (at resonance), P starts to saturate (cf. fig. 1.3) and reaches 
1/4717 


O h2 l 2 3 


Fig. 1.3. Pm and yo/z" as a function of h?. 
Pm ieee ly ae 


asymptotically the value 
N(hw)?S-(S+1) 


(1.16) Po = 3hT 


At 4° K and H=3000 6 typical values for the relaxation constants 
are tT ~ 10-98. and 1 ~ 1073s. requiring a r.f. field of about 0.1 @ to 
reach half the expected saturation. 


[1.3] Microscopic theory 


We shall now approach the subject from a microscopic viewpoint. 
Suppose we have an assembly of spins of which the spin Hamiltonian 
(1.04) is known, and that the eigenvalues and eigenfunctions have been 
calculated. 

When a radio frequency field hz sin wt is applied as a perturbation 


(1.17) H .=yhhz(o)Sz sin wt, 


then the transition probability induced by this perturbation will be 
proportional to the square of the modulus of the matrix elements 
<i|Sz|7>, taken between the eigenstates |i> and |j> of the spin Hamiltonian. 
i —E 
; j 
to spin interactions and lifetime limiting processes. The shape of the 
transition will be characterized by the normalized shape function g(@— wij) 
where 


will have a certain width due 


The transition frequencies wij = 


+00 
dea 
(1.18) | g(o — oy) 2 = 1. 
—oo 


If the transition frequencies have a separation much larger than 
their width, the induced transition rate Wi; can be written as 


(1.19) Wig = /ay?h?(w) g(o — 49) <2|Sz|)>|? 
with 
(1.20) Wiu= Wi. 


These results follow from normal time dependent perturbation theory 
(cf. [10]). 

If the population of the ith level is called n;, then the power absorbed 
by the spin system due to transitions from i to j will be 


(1.21) Pm=(nj—ni)hoiy Wi. 


The populations will change, partly due to transitions induced by the 
v.f. field and partly due to transitions induced by the lattice vibrations. 
(Spontaneous emission can be neglected.) 

We shall now consider a two-level system and call the lattice induced 
transition rate Uys, level 1 being the lowest in energy. 
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Multilevel systems will be considered later as they lead to rather large 
formulae, masking the effects we want to discuss here. 
For the rate of change of the population difference we find 
d n1—Nneo yr , 
(1.22) —— = 2(—nU 2 + 22U21) — 2(n1 — na) Wyo. 
Because of the fact that the energy exchange between the spins and 
the lattice, which is in the thermal equilibrium at a temperature 7, 
takes place in quanta fiw, we have 
(1.23) Ure _ ot(horw/kT) _ Me 
Uar ny 
where 7) and 7. are the equilibrium populations, having their normal 
Boltzmann ratio. 
The exponential may be expanded as in most cases hai, < kT’, and 
(1.22) can be rewritten using the relation nj+n2=1)+ M2 as: 
d(n1 —N2) 


(1.24) a [(m1 — n2) — (M19 — N29) ] 2012 — (m1 — 22) 2Wite 


where 
ie. _ Ux = Uar 


The steady state solution of (1.24) is 


N11, —N2, 


(1.25) n—n2, = ea 


If N is the total number of spins we find with (1.19) 


N Awe 1 
Ny — N29 = IS+1 kT and ioe = et 
howWie _ 
\ Pn neo) Wain = 
2 © 
(1.26) (hier)? 1 /ay2h?() gl —eo12) {C1 Sz12>/2 
~ 28-1 kL 141 /ey?h?(w) g(w@ —aie)t1|<1|S'2|2>|? 


[1.4] Spin-spin equilibrium and the resonance line shape 


[1.4.1] Spin-spin equilibrium. There have been many theories 
and experiments on the internal equilibrium in spin systems. The main 
body of the experimental work has been carried out by GorTER e.a. [11], 
using their relaxation methods. The spin system has many degrees of 
freedom, more or less coupled through the spin-spin interactions. If such 
a system is in internal equilibrium one can assign a temperature 7's to 
it and by using its density matrix 

exp —#s/kT's 
OS™ Yexp —H s/kT's 
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the equilibrium properties can be described. In most cases one does not 
have to take the total density matrix of spin system plus lattice because 
the spin and lattice degrees of freedom are independent, apart from a 
small term causing spin lattice relaxation. To obtain gs one can sum 
the total density matrix over all lattice states. 

That the spin system is in internal equilibrium does not merely mean 
that the populations of the energy levels are in accordance with the 
Boltzmann factors. It means also that the non-diagonal elements of the 
density matrix are zero. Or, macroscopically, the transverse magnetization 
is zero, 

Let us consider what happens macroscopically in the Bloch picture 
when at a certain time the magnetization M is not along the magnetic 
field (cf. [9] and fig. 1.4). Because of the coherence disturbing effects the 


Fig. 1.4. Time behaviour of M according to Brioc#. 


transverse magnetization decreases according to (1.09) with a time 
constant tz. During this process the internal energy M-H remains con- 
stant, as in the Bloch picture a change in the spin-spin interactions is 
neglected. After a time of about t2, M will reach the value MM, and the 
spin system will be in internal equilibrium at a temperature 7's which 
corresponds to its internal energy and the size of M,. 

According to Curie’s law (1.08) we have 


(1.27) M; = = where a = ae 

If in paramagnetic resonance the r.f. field is much smaller than the 
coherence disturbing internal fields, it is still legitimate to describe the 
situation by means of a spin temperature. To describe the saturation 
behaviour one has to put in (1.26): 


(1.28) Ts=T (14 1oy*h?g(w— iy) t1|<i|Sz|j>|?]. 


At r.f. fields higher than the internal fields, or for times shorter than T2 
after a shock excitation of the spin system, it is not possible to speak 
of a spin temperature. Phase coherence effects then play an important 
role and the state of the spin system cannot be described in terms of 
phaseless constants. 


The studies of GorTER e.a. in this field have been mainly concerned 
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with the spin absorption at low frequencies which is closely connected 
with resonance line shapes, although no completely satisfactory analysis 
has yet been presented to relate them [12] [13]. Especially in the case 
of multilevel spin systems, as will be discussed in [1.6], and systems 
with exchange the theory is far from complete. 

There have been many discussions as to whether it is allowed to 
split the spin system into a Zeeman part and an interaction part each 
with its own internal energy and temperature. This procedure has been 
followed by BLOEMBERGEN e.a. [14], Yoxora [15] and ABRaGame.a. [16]. 
The difficulty is that the degrees of freedom are not independent. In the 
case of large exchange interaction they may be considered to be commuting, 
which is not sufficient to justify their assumptions. 

According to a theory of Kroniag and BouwKamp [17] the rate of 
spin transitions due to the spin interactions in a two-level system, is 


(1.29) Us = Wd Ee 2/24" 


where ho is the energy required for the spin transition and wg a frequency 
characterizing the dipole-dipole interactions, which will be introduced 
later (cf. (1.35)). One may regard the exponential as the probability that 
the internal field just compensates the applied field, or as the amplitude 
of the Fourier component at frequency of a gaussian random fluctuating 
internal field. This probability readily becomes a very small quantity 
as mw is increased, one of the reasons why this theory was originally 
rejected by Brorr [18]. 


[1.4.2] The resonance line shape. Comparing the results of the 
quantum theory and the classical picture, it is found that (1.26) corresponds 
with (1.15) when 


22 


NO OL) = To —arere 


This is a Lorentz line shape with a half width at half power of 1/t2. It 
is not surprising that the classical description leads to a Lorentz line 
shape, as it was introduced by assuming (1.09). This is one of the main 
drawbacks of the Bloch picture as experimentally one finds mostly an 
almost gaussian line shape in the absence of exchange interactions. 

The shape of a paramagnetic resonance line will partly be determined 
by the lifetime broadening due to spin lattice relaxation. That contribution 
will have a Lorentz shape with a half width 1/1. There are many salts 
in which this is the main contribution to the line width at the higher 
temperatures e.g. TiCs alum [19], FeNH, tutton salt etc. 

This situation arises also in salts containing ions of the rare earths 
group, but for the iron group we have for the ions with an odd number 
of spins in general 1; <2, making the lifetime broadening negligible. 

To calculate the contribution of the spin-spin interactions to the line 
width two widely different methods have been used. VAN VLECK applied 
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the method of moments to the case of identical spins [20]. He calculated 
J (@—or)4g(@—a@z)dw and f(@—or)29(@—a@z)dw 


by a diagonal sum and showed that #,, has no influence on the second 
moment of the absorption line, because of the commutation properties 
of #,., Hz and #,,. There are two contributions to the fourth moment, 
however, the first one due to dipole-dipole coupling alone, the second 
due to the rate of change of this coupling which, as we shall see later, 
is determined by *#,,. 

ANDERSON [21] and Kugo and Tomrra [22] have developed a mathe- 
matical model which allows one to calculate the actual line shape. They 
assume the system to absorb a single frequency, which varies in a random 
way over a certain range determined by the dipole-dipole interactions, 
but at a rate controlled by the motions or exchange coupling. This concept 
was first proposed by GorTER and VAN VLECK [23]. 

We will first split up the dipolar interaction (1.02) into terms depending 
on their influence on MW, 


(1 3 - ) H a = > AS, Sz; ' BUSTS; - S;_S7*) = C(S7 Sz; + SzSj*) + 
.30 i>j 
+ DS; Sz; + SzS;) } EB SrSF + FS;S; 


where St=S,-+ iSy and A, B,C, D, E and F are dependent on the polar 
coordinates of ion 7 with respect to ion j. 

ABRAGAM Called these terms stop-stop, flip-flop, flip-stop ete., according 
to their effect on M,. 

For large static fields the terms A and B will broaden the absorption 
line at the Larmor frequency wz, while the terms C, D, # and F introduce 
weak satellite lines at frequencies around 0, 2m, and 3ez (ef. fig. 1.5) [24]. 
The absorption near zero frequency has been the subject of many experi- 
mental studies by GorrmR e.a. [11] [25] [26]. 


* 
xtenyfarl RY 


——— 
/ 
QO 7, 


Fig. 1.5. 4”(w)/m as a function of w/ar. 


Ma < Wy < WL eae wl <= Wa. 


[1.4.3] The influence of exchange coupling. A brief out- 
line of Anderson’s theory [21] [27] will now be given. 
We take as the Hamiltonian 


(1,31) MS GEE Nets ne > Fix S7-Se+gBh > Say. 
j>k i 
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The absorption spectrum as measured by the last perturbing term is 
given by the Fourier integral of the time dependent operator S(t) 


(1.32) I(w) = Trace |f Se(t) edt? 
where 
(1.33) ibe — (0, Sa), 


Since #,, commutes with #z and with );S,,, the exchange interaction 
causes no broadening as it does not directly change the energy levels nor 
interrupt the radiation, but we do have 


(1.34) th Ha=(H24+H q, Hal ~ (Ho, Ha: 


The omission of #z in the commutator expresses the assumption that 
Hq is so small as to have no appreciable matrix elements between the 
eigenstates of #z and also that we are only interested in the line near 
the Larmor frequency. 

This assumption is the much discussed truncation and only valid when 


KH p@§ KHz and Hay KH 2. 


We thus retain only the secular terms A and B of #q which results in 
a frequency shift of the resonance. This shift has a time variation due to 
(1.34) at a rate of the order of #%,,/h, to the effect that Hq is partly 
averaged out during the time that resonance is observed. 

The assumption now is that this time variation is a gaussian random one. 

The results of calculations on this basis are that the second moment 
has for a simple cubic lattice the value 


(1.35a) Og= <l@ == wr|?),,= 117433 yn2s : (S+ 1) (Ai4 + Agt+ A34 = 0.187) cd-6 


and the fourth moment for H along the (100) axis. 
0.021 Tt (Wex\2 
(1.35b) <|o—oz|*>ey = Bwa4| 0.742 +e-1{0.098 | a 5(2=) ] 


where d is the lattice constant, c the magnetic concentration and Aj,2,3 
the cosines fixing the direction of H with respect to the crystal axes. 
The half width at half power of the resonance line becomes 


Aw = ee for waX x. Kx 
(1.36) Dex 
Aw =oq for Mx <a 


with 
8,48 /J\2 
(1837) = == (7) S(S41). 


For the case that w,, > oz the truncation performed in (1.34) is no 
longer permitted as the energies for transitions induced by the non- 


23 Series B 
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secular terms of #q are now available from the interactions. In the 
calculation of #q retaining all the terms of (1.30) multiplies (1.35a) by 


the factor = 


We get 


(1.38) Aw =~—— for wo, > or. 


[1.4.4] Inhomogeneously broadened lines. It may be that, 
upon diluting a substance magnetically, #q becomes smaller than 
variations from ion to ion in the zero field splitting or static local fields, 
due to ions of another kind or field inhomogeneities. For the so-called 
inhomogeneously broadened lines then occurring one has to modify the 
absorbed power (1.26). For g(@—«jy) the single ion absorption shape is 
taken, but as the ions are spread over a wider line, one has to multiply 
the whole expression (1.26) by DS dah 

“  g(@—orz) 
g*(@— ao) to the single ion shape, and integrate over the distribution of 
the Larmor frequencies of the different ions. If g*(@—wo) only varies 
slowly over the region where g(@—qz) is appreciable, then the magnetic 
absorption becomes: 


, the ratio of the observed shape 


(1.39) Pit A el ee eee oe 


2S+1 kT V1+/ey*h?g(0)t1<1/Sz2]2)2 
where g is assumed to have a Lorentzian shape [28] [29]. 


(To be continued ) 
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[1.5] On the theory of the relaxation rates 


The energy transfer from a spin system to its environment has been 
the subject of many controversies, and the different mechanisms that 
have been proposed will be reviewed in this chapter. We shall limit our- 
selves to the case of non-conducting ionic crystals. 


[1.5.1] Spontaneous emission. An excited spin can return to 
its ground state by spontaneous emission, in which case it gives a quantum 
hw to the radiation field. For the case of a single excited spin in free 
space the Einstein coefficient A for spontaneous emission is: 

3 
(1.40) A= “ | Mmn|? 
where /imn is the matrix element of the component yhS, of the magnetic 
moment perpendicular to the static magnetic field H. At frequencies w/2x 
around 1019 Hz the value of A is about 10-! giving a lifetime of 
1012s. or 105 years. 

In resonance experiments, however, the conditions for the calculation 
of A are not fulfilled, and this rate will be enhanced for two reasons. 

The first is that the spins are not in empty space but in an electro- 
magnetic resonator with a quality factor Q@ and a volume V¢. PURCELL [30] 
has pointed out that in this case the energy density of the field per unit 
frequency range and unit volume is increased by a factor B=/3Q/82?V¢ 
making the rate for spontaneous emission 


4nQ F 
(1.41) Wea ae Ve | mn|?. 


In our experiments 6 ranged from 50 to 500. 

The second reason is that we are not dealing with a single spin, but 
a great number of spins, constricted within a space with dimensions 
smaller than the velocity of light c divided by the frequency width Aw/2z 
of the transition band. In such aggregates of spins coherence effects 
become important and states exist with radiation rates far greater than 
normal [31] [32]. However, the factor by which the radiation is increased 
is for paramagnetic materials in saturation experiments of no importance. 
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[1.5.2] Lattice induced transitions. For the lattice vibrations 
there are two ways to influence the magnetic moment. One way is through 
modulation of the spin-spin interactions, the other through modulation 
of the electric crystal field near a magnetic ion. 

In spin lattice relaxation, one spin quantum is used to excite lattice 
states, and as these states have a continuous spectrum, one can apply 
time dependent perturbation theory; the calculations differ from spin- 
spin relaxation ones, in which one cannot describe the relaxation mecha- 
nism in terms of a coupling between one spin and the spin system as a 
whole. 

WALLER [33] has pointed out that two essentially different processes 
are possible by which the lattice vibrations can take or give up the energy 
for the spin transition. When the spin jumps from level 7 to level 7 the 
energy 


(1.42) E,— Ej =howi 


is in the direct process compensated by the emission or absorption of one 
elastic quantum ’|@,j| and is induced by the part of the spin Hamiltonian 
linear in the displacement vectors uz, the term #1, in (1.07). 

In indirect processes two lattice quanta play a rdéle. One ha; is absorbed 
while another hw, is emitted so that 


(1.43) h\ w4j| =hox —hor. 


This is called a quasi Raman process and has its origin in the quadratic 
terms #2,, in (1.07). 

The lattice will be considered to behave according to the Debye theory. 
We are interested in the part of the frequency spectrum below 1/200th of 
the Debye temperature when dealing with the direct process. So this 
assumption is not a real restriction. For the indirect processes the whole 
frequency spectrum plays a réle, and so the optical branch as discussed 
by BLacKMANN [34] may have an important influence at the higher 
temperatures. 

In the Debye theory, according to the procedure of SoMMERFELD and 
Berue [35], the normal modes of the lattice are quantized. The lattice 
states are described by a set of quantum numbers p,; for each of the kth 
modes of the elastic vibrations. The energy of the mode is given by 


(1.44) Ex= (pe+3/2)hor. 


The normal coordinates gz will have matrix elements 


(1.45) <pe|ge| pe + 1>= | a )" 


where M is the total mass of the crystal. As VAN VLECK remarked [36] 
this is a weak point in the theory as it assumes that the amplitude of 
vibration is the same for all atoms regardless of their mass, although in 
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first approximation one can neglect the particle structure for long phonon 
wave lengths. 
The number of oscillator modes in a frequency interval between w and 
w-+dw for the longitudinal waves with velocity 1 is given by 


w2Vda 


270,83 


(1.46) oi(w)do = 


and for the transverse waves with velocity v; by 


w2Vdw 
272v,3 


(1.47) o:(@)dw = 


where V is the volume of the crystal. These formulae are valid for 
@ <We, the cut-off or Debye frequency. 

If N is the number of atoms in the crystal, there are N longitudinal 
and 2N transverse modes of vibration, so 


o @ 


f oi(@)dw=N and eon 
giving j 

(1.48) Oq = (622 5)" and We; = (6x2e 7) 

We can write the displacement vectors u; as a linear combination of 
the normal modes qx of the lattice, 

u;= 2 An Gx Sin wx(t —7;/Va+ Vx) 


where vg=v for the longitudinal modes and equals x; for the transverse 
ones. 

The relative displacement between the atoms 7 and 7 for waves propa- 
gating along ry, if the wave length A, > |ri|, is 


(1.49) Uy =U; — Uy = Tj 2 Ok — qx COS ox(t—r/va+ x). 
k 


For the direct process we are interested in matrix elements of the linear 
part of (1.07) corresponding to the creation or destruction of one phonon: 


(1.50) <p, t|#1(qx; Si, Sm)|px+ 1,7> 


with the condition (1.42) due to conservation of energy. 
The rate Uy, at which transitions are induced in the spin system by 
this direct process will be 


(1.51) Uijg = F olox)|(2#1) 2dr oe f (heox)*|(1'>|2g(on — oy) (Pe + 1) gyda 


where the exponent « gives the frequency dependence of @(wx)|<# ba tar 
and <#1'» contains the matrix elements of the spin operators in <i|41|)>. 
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The integrant of the second integral, which extends over all possible 
lattice vibrations, contains the distribution function of the spin transition 
centered at wij. 

If the lattice is in equilibrium at a temperature 7’, then the average 
number of phonons at a frequency is 


oe L 
9 Nr —— 
Leen) Pk av ~ Tok 


Combining (1.51) and (1.52) one obtains for wij > wa 
(1.53) Uijg & wf-1-T for hay <kT. 


thus at very low temperatures, the zero point vibrations determine the 
relaxation rates. 

For the indirect process we are interested in the matrix elements of #2 
in (1.07), quadratic in the normal coordinates: 


(1.54) Pk Pts | H2(Gu, qu, Sm, Sn)| pe +1, pi— 1,7 


with (1.43) as a condition. The transition rate due to this indirect process 
will be 


@ 


(1.55) ain = re 0(@x) o(@2) |X H#2>|? g(w@xE — wr — wiz)daox dwt. 
0 


The main contribution to this integral will come from phonons around 
kT or hwe, depending on which is smaller. One has to go to a higher order 
perturbation calculation but the decrease in the matrix element is more 
than offset by the large increase in the number of lattice vibrations that 
can take part in the process at higher temperatures. 

We substitute m,—@:=q@ into (1.55) and find: 


@, @, 
(1.56) Uitina © Jf og g(@ — wiz) |(H2">|? (e+ Day Pe >ay Un deo 
0 0 


where (#2'> contains the matrix elements of the spin operators in (#2). 
With (1.52) one finds for (1.56) 


(1.57) Vijina & T? when kT S hare 
(1.58) Uisna & T*+1 when kT < have. 


The case wij> me can be neglected. 


[1.5.3] The energy transfer by means of modulation of 
the spin-spin interactions has first been considered by WALLER [33] 
for systems without exchange interactions. One may include the exchange 
interactions by considering the magneto-elastic interactions as a pertur- 


bation which causes spin transitions because of lattice strains (magneto- 
striction), 
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The dipolar coupling #@ (1.02) depends on the radius vector ri; between 
ion ? and j and will therefore be affected by the lattice vibrations. 

In (1.07) the linear part becomes proportional to yh lriyz4 and the 
quadratic part to y2h2/r;;5. The non-secular terms of these perturbations 
come from the corresponding ones in #q in (1.30), and are derivatives 
of the terms C, D, E and F, which can cause transitions in the spin system. 
WALLER arrived at values for the relaxation rate of about 10-1 s— at 
liquid air temperatures and 10-4. H?s-1 in the liquid helium range. 
As the observed values for a typical salt such as CrK alum are 10% s-! and 
102s resp. this mechanism is too weak to explain the observed values. 
When there are inequivalent ions in the unit cell, as is the case for most 
salts, the terms A and B will also contribute to the relaxation. This, 
however, does not much reduce the large discrepancy. 

The exchange integral J is a sensitive function of the radius vector rjg. 
If ions 7 and k are equivalent, however, the exchange interaction between 
them commutes with the Zeeman energy and #,, will not be able to 
change the magnetization. 

In [1.4.3] and the subsequent discussion we have seen that the exchange 
produces a random frequency modulation of #g, and one could therefore 
imagine a perturbation #z—> Hq > #,, — lattice. I have made an 
estimate of the relaxation probability due to this process, assuming that 
the exchange constant J; in (1.03) is an exponentially decreasing function 
of the distance riz between 7 and j, and that ,, < wx. 

The amplitude of the Fourier component at frequency wz in #q, due 
to the modulation of #,, by the lattice vibrations, can be calculated by 
using (1.34), (1.49) and the above-mentioned assumption. The effect of 
this Fourier component on the spin operators is found by using (1.33). 
For the direct effect it is found that Uij,« 3,7 where ,, is defined by 
(1.37). For w,, ~ 3.1010s-1 a rough estimate makes Uy, ~ 10-> s+ at 
1° K. As this rate is quite slow, this mechanism seems not to be the most 
important one. 

BLOEMBERGEN and WANG [14] indicated that another situation arises, 
however, when the exchange frequency (cf. 1.37) is much larger than 
the Larmor frequency. For this case the 10/3 effect comes into play. 
The energy for a spin transition is now available in the dipole term 
because of the strong modulation by #%,,. An analogous situation arises 
in the case of nuclear resonance, where the dipole term is modulated by 
the random velocities of the nuclei. Up to a cut-off frequency @ one 
finds a white spectrum there and 1 becomes about equal to t2 if @ > wr. 

When o,, > wr the energy exchange between “exchange system” and 
lattice is at a high rate, as it can take place in large quanta hw,,. The 
bottleneck for the energy transfer will presumably be between the “Zee- 
man” and “exchange” system. With respect to the random frequency 
modulation picture cf. [1.4.3] and the above-mentioned analogy, BLoEM- 


BERGEN e.a. expect T] * Te. 
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An example of this case is the free radical diphenyl pycril hydrazyl. 
BLOEMBERGEN and Wang report a value of t of 6.4 10-8s. at 300° K 
and 77° K with t2 ~ 6X10-8s, while Meyer [37] found 11 ~ 10~’s. at 
4.2°K, VERSTELLE (to be published) has recently studied this effect in 
more detail. 

As this type of relaxation is determined by the spin-spin interactions 
and independent of the lattice vibrations one would expect U not to 


depend on temperature. 


[1.5.4] The influence of spin orbit interactions. Xelaxation 
probabilities due to the modulation of the crystal electric fields are 
essentially single ion processes. HerrLeR and TELLER [38] considered the 
case of substances with a spin splitting at zero magnetic field. The lattice 
vibrations modulate the value of D in the spin Hamiltonian (1.04) and 
induce transitions in the spin system accordingly. They considered only 
the direct process. Frmrz [39] extended these calculations to the case of 
indirect processes. Their results are similar to those of WALLER’s when 
gpH is replaced by 2D and they find for the direct process 


Uy ew 3.3 x 106 wi2kT hag. 


Kronie [40] was the first to recognize the important possibility for heat 
transfer via the orbit lattice coupling. VAN VLECK [36] independently 
carried out similar, more detailed calculations for the special cases of 
Cr and Ti alum. 

In the alums the paramagnetic ions X are surrounded by a near 
octahedron of watermolecules Y. 

Van VLECK [41] considers the influence of the normal modes of vibration 
Qi of the cluster XY, and their effect on the orbital momentum L of X. 
The modes Q; can again be developed into linear combinations of the 
lattice modes qj. 

The crystal field of the cluster XYg around the center ion X will 
generally contain only terms quadratic in the Q’s. Our case differs, how- 
ever, from the general situation in two ways. If an orbital state of the 
central ion is degenerate, the Jahn-Teller theorem [42] states that the 
octahedron will become distorted, so as to lift this degeneracy. The atoms 
around the cluster will partly oppose this distortion. Through this Jahn— 
Teller effect terms linear in the Q’s with even symmetry will appear in 
the field. 

Another way to obtain linear terms is by the so-called indirect effect 
of atoms not belonging to the cluster X Y¢. In an early analysis Van VLECK 
[43] comes to the conclusion that the direct effect of distant atoms is 
not very important, but the influence of them on the octahedron of water- 
molecules is quite large and will bring linear terms in the field acting on X, 
More recently many paramagnetic resonance experiments have shown 
that this indirect effect can be large. The constants in the spin Hamiltonian 
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of a paramagnetic ion are found to depend on temperature, magnetic 
dilution, pressure and the size of the other ions present. It may be that 
this indirect effect of distant atoms is coupled with the Jahn—Teller effect. 
Only the even normal modes produce linear terms in the crystal field. 
They are (cf. [41]) Qi (total symmetric vibration, representation I, ), 
Q2 and Q3 (tetragonal distortion, P'3 a Qa, Vs and Qe (trig. dist. I’s,) and 
Mis, Goo and Qe: (rotation, irae): 

Considering (1.04) and (1.05), there are two different ways for the 
lattice vibrations to influence the g-factor: 


A. The adiabatic terms. 


These are present when considering the change in the orbital wave 


s |n> by the Q’s only. The perturbation calculation can be carried 
out by keeping S and Q as operators, but as (H,—o) in (1.05) retains 
the same value, the second order perturbation becomes zero because of 
the Hermitian property of any element and the pure imaginary values 
of the matrix elements of # so, while those of #oz are real. 

VAN VLECK [36] showed that one has to go to a third order perturbation 
to get non-vanishing terms. A good example of these terms being pre- 
dominant is Cr+++ in CrK alum. In equilibrium the water octahedron is 
slightly distorted along a body diagonal of the cubic alum crystal, thus 
producing a crystal field of trigonal symmetry with this diagonal as a 
symmetry axis. By taking the 2, y, z-axes along the principle cubic axes, 
Van VueEcK finds for the interaction between S$ and Q 


H s1= €1(Q3( 2822 Sz? —Sy2) + V3 Qa(Sy2—Se2)] + 


(1.59) 
+ 2 [Qa(S2Sy “i SySz) ahr Os( az ar Sz) ail Q6(SySz te S2Sy) | 


where ¢, and eg are constants proportional to /2/(Hy—Ho)?. 

After expressing the normal coordinates of the cluster in the normal 
modes of the lattice and multiplying the square of the modulus of the 
matrix elements by the phonon density, the relaxation rates become: 


(1.60) Uigg=1/h?<0(w1y)| 4 sx (tpn, (Pr + 1)) | >ay- 


With (1.51), (1.52) and (1.53) one has Uy«Tij?. At low fields and at 
TJ =1.4 °K it is found that U » 50s71. 

To obtain the power which the spin system is able to transfer towards 
the lattice from the values of Uy, one has to take a suitable average over 
all energy levels. 

We will postpone this calculation to [1.7] where a more general discussion 
of the power transfer will be given. 

The values of Ujj,,, calculated with (1.55) are independent of wi. 
The value of the exponent « for the adiabatic processes equals 6, made 
up of an wt dependence of the number of lattice vibrations taking part 


in the process, and an w? factor due to (1.49). 
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One of the underlying assumptions for calculating Uy,,, is that the 
phonon wave length is large compared to the dimensions of the water 
octahedron, a condition which is only fulfilled at temperatures well 
below the Debye temperature, which is about 80° K for the alums. 


B. The non-adiabatic terms. 

If one takes into account the contribution of the energy of the phonon 
or spin jump hay to the energy difference in the denominator of (1.05), 
the second order perturbation of the lattice vibrations on the spins does 
not become zero, and can be relatively large in cases where the energy 
difference of the orbital levels between which L has non-vanishing 
matrix elements, is small so that one cannot neglect ha compared to 
E,— Lo. 

An example of this case is TiCs alum where this splitting is about 
400 cm-! which is small compared to the corresponding value of about 
23000 em-! for Cr+++ in CrK alum. 

When developing the energy denominator, in which because of energy 
conservation the phonon frequency has to be wi , one finds for the matrix 
element of #s, by perturbation with the orbital levels n 


<0, 1, pel AH sx|0, 7, pe + 1> = 
1.61 ae SO ee es | 
sti ape 2 i. Eo (0,1|#so|n,7> <n, pe| Hor] 0, pe + 1>. 


The relaxation rate will have the same form as in (1.59), but due to the 
appearance of wy in (1.60), one finds Ug « wijt*T for hay <kT. 

For the indirect processes the exponent « in (1.56) will be 8, which 
is 2 higher than for the adiabatic terms for the same reason why «wi; 
appeared in (1.61). 

For Cut++ in CuK tutton salt calculations have been made by NAGAaoKa 
[44] on the same basis as those of VAN VuiEcK. The five orbital levels 
in the lowest multiplet are, in a cubic crystal field, split into a lower lying 
doublet and a triplet with a separation A of 15500 cem-, 

The rhombic component in the electric field splits the doublet with a 
separation of say A’. This last splitting has no influence on the g-factor 
in first order, but does influence the relaxation rates. The adiabatic and 
non-adiabatic terms seem to be of about equal importance. NAGAOKA 
expects a highly anisotropic Uq« H4 depending on (A’/A)?. This rate is 
expected to become extremely small when H is parallel to one of the 
axes of the crystal field. 


[1.5.5] The Temperley effect. It has been known for a long 
time that the discrepancies between these theories and experiments are 
quite serious in the liquid helium range. The relaxation rates in that 
range for the concentrated salts are usually found to be a decreasing 
instead of an increasing function of wy while the temperature dependence 
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is stronger than expected for the direct process, but not strong enough 
for the indirect ones. So for undiluted CrK-alum it is found that U «78 
and for undiluted FeNH,-alum U «75. The diluted salts seem to behave 
more according to VAN VLECK’s theory. 

TEMPERLEY [45] attempted to explain these discrepancies by supposing 
that by spin-spin interactions a number of spins jump at the same time, 
creating a phonon with the total energy. This effect could explain the 
field dependence of the relaxation probabilities. To explain also the 
temperature dependence one has to assume that a large number of spins 
jump together making the excited phonon frequency w > k7/h. 

Van VuEcK [36] doubted whether this effect is large enough. We 
shall discuss it later in more detail. 


[1.6] Equilibrium and saturation of multilevel spin systems 

[1.6.1] The attainment of equilibrium in spin systems with 
S>1/, can be different from the case of S=1/2 in this respect, that in 
addition to the processes considered there, the separations of the energy 
levels may be such as to allow flip-flip or flop-flop processes requiring no 
energy. These processes become possible because of energy splittings at 
zero magnetic field. 

GoRTER, DREWES and VERSTELLE [12] used this picture for the 
description of their experimental results, but did not work it out so that 
it explains the experimental results quantitatively. 

BLOEMBERGEN e.a. [46] considered the possibility that a multiple 
effect of the spin-spin interactions is important. The idea is that at points 
where the resonance lines of different spin transitions overlap, flip-flop 
processes (« B2in (1.30)) between neighbouring ions occur. They call the 
rate at which this takes place the cross relaxation rate Ur. 

Jumps in which more than two spins take part, may sometimes increase 
Uer considerably, and for LiF it is possible to explain the short times 
required to establish spin equilibrium [47] by considering triple spin 
processes. 

For the overlap they take the product of the shape functions of 
the two resonance lines. If these lines are far apart compared to the line 
width, this overlap will be small, reducing the cross relaxation 
rate considerably. 

When saturating a single resonance line of the absorption spectrum 
Ucr becomes important when comparable with the spin lattice relaxation 
rate. Calculations for similar processes have been made by ANDERSON 
[21] when he considered the effect of the exchange on the fine structure 
of resonance lines. The jump rates he considered were of the order of 
wq. At this moment we are interested in much smaller rates, of the order 
of 10-% wa. 

The resonance transitions can be considered as weakly coupled oscil- 
lators. In the case considered by ANDERSON the coupling is strong 
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and one therefore gets an appreciable distortion of the initial resonance 
lines. In the case considered here, however, the distortion will be too small 
to be observable. The exact values of U;, will be difficult to calculate, 
but there is no doubt that in concentrated salts the interactions are 
usually so large as to produce saturation of the total spin system when 
saturating a single resonance line with a radio frequency field. 


[1.6.2] No cross relaxations. We will consider first the case of 
low magnetic concentrations without cross relaxation. At times much 
longer than the spin-spin relaxation time t2 after applying a r.f. magnetic 
field to a sample in a static magnetic field, the transverse magnetization 
will be damped out, and the spin system can be described by the population 
densities n; of levels 7. There are 2S+1 levels and we have always 


(1.62) Y m=N 


the total number of ions present. We shall discuss the rates of change 
of nj and the influence on n; of radio frequency fields at the various 
resonances. This can be done by setting up equations of the type (1.22). 
The transition rates from level 7 to level 7 due to interaction with the 
surroundings are called Uj; and those due to r.f. fields Wy; the sum being 
Vi= Ui; + Wij. 
Because of (1.23) we have 


(1.63) Uyg= Uy ehoy/kT 


with haiy=H;—H; and also we have Wiy= Wy. 
We will now use a treatment analogous to that of LLoyp and Paxz [48]. 


(1.64) a LY’ = iVig try Vii (1,7=1...(28+1)) 
] 

where the prime means that the sum is taken for all values of 747. Of 
the 2S+4 1 equations (1.64) there are 2S independent ones, from which 
one variable can be eliminated with (1.62). We are left with 2S simul- 
taneous differential equations with 2S independent variables. 

With a trial solution m—Ni,=Are~* for the equations for the deviation 
from equilibrium one will therefore in general find 2S time constants. 
If nig is the equilibrium solution we have in general 


28 
(1.65) Ni = Nig + > Anew! (= 1,2,...28+1) 
k=1 
where nip, Ay, and A, are functions of V vg: 
If after some time several terms in (1.65) are dominant, it will not 


be possible to describe the population differences with a single relaxation 
time. 
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Experimentally one can measure the population differences between 
two levels by measuring the resonance signal (1.21). 

The steady state solution of (1.64) is obtained by putting dnj/dt =0. Being 
interested in the behaviour of the n,;’s as a function of the induced transition 
probability between two specific levels, for which we can take levels 1 
and 2 without loss of generality, we shall write out the equations for this 
case, substituting n1—n2=Ale 


—Aje > (Via3+ Uy+ W 12) + n2(Ue1— Ui2— > Vij) +... 


j+1,2 j+1,2 


a 
3 
cS 
—_ 
: 
cS 
S) 
=i 


» + Ness Vess1,1=0 


Aj2(U12 + Wiz) +me(U12—Un— > Voj)+... 


2 
+NpVpq +... +Nes41V e541, 7=0 


Aj2Vi,2s+1 +ma(Vi,ess1+ Vo,es41) +... 
+npVp2su+...+nest > —Vesis,7=0. 
G+2S+1 
Wi, appears in the first column in the first two equations. One of them 
can be eliminated by replacing for instance the first equation by (1.62) 
and we then have 2S+1 independent equations with 2S-+1 unknowns. 
With the well-known rule of Cramesr for solving a set of inhomogeneous 
linear equations Aj2. and ng can be found. Of the determinant D of the 
coefficients, with Wy. put zero, we call Cpg the cofactor of the pth element 
in the gth row. The solution is 
NCu 


(1.67) Ai2 = DiWia 


NC + Wi2Cx1, 12 


(1.68) DE aa ag 


Cia,12 is the cofactor of the term Ui2+ Wiz in Cr. 

Analogous to (1.25) and (1.26) an average relaxation time can be 
defined characterizing the saturation of Aj2 
= OT 2? os ti Vir 
1-69) T12 = a a O12 
For the difference between the populations of level 7 and k one can write 
with (1.68) 


Oj1 —Cr + Wie(C jr, 12 —Cr, 12) 
(1.70) Ae = SDA WvaCam a 


The powers absorbed at the various frequencies in the case of well- 
separated resonances is found from 


(an) Py =Ayhoy Wi. 
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Relations between the various cofactors can be obtained by realizing 


the in equilibrium 
a N AWK 
(1.72) (Ajx) w=0 = 3841 kT’ 


The cases pertaining to the experimental conditions will be calculated 
in the appropriate places. 


[1.7] The power transfer between spin and lattice systems 


The relaxation rates as discussed in [1.5] enable us to calculate the 
power transfer between spin and lattice systems. For the case of cross 
relaxation rates faster than the spin lattice relaxation rate, the spin 
system will come into equilibrium at a temperature 7's in a time 1/Uer, 
short compared to 11. The spin lattice relaxation will then take place 
according to a simple exponential decay e~“", and is due to the power 
transfer of all transitions. To express the power transfer between spins and 
lattice we will not use the populations n; of the different levels. This can be 
done for the case where the spectrum consists of well-separated lines, 
but if these are overlapping ones, individual level populations lose their 
meaning. 

For the power absorbed from a r.f. magnetic field by the magnetic 
material, we have (cf. (1.26)) 


(1.73) Pm=QW/T's 


where W is the induced transition rate Wj4;, summed over all possible 
Larmor frequencies present in the material. The distribution of the spins 
over their energy levels is described by the density matrix. 

cp —E;,/kT. 
1.74 Ey) = es —— 
( ) os( i) di exp—Hi/kT's 


where H; are the eigenvalues of the spin Hamiltonian #%,, with inter- 
actions, 


[1.7.1] Direet processes. For the lattice induced transition rates 
due to the direct processes, we take the integrant of (1.51), without the 
shape factor. Integration over the shape of the transition as was performed 
there, is not necessary as we include here the spin-spin interactions in 
the spin Hamiltonian. 

The energy transferred from the spin system in temperature equilibrium 
to the lattice, by the direct process will be: 


(1.75) Pe = LY hois|Uijaos(Hi) — Ujiaos( #5). 


Due to conservation of energy we have m=; and with (1.51) formula 
(1.75) becomes 


(1.76) Pe, = 2 (heoiy)**1| <i] 1’ (S)|7>|? {es(Hs)(py + 1) —05(E;) pu}. 


ij 
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The factor in brackets contains the part dependent on spin and lattice 
temperature T,. 

For T's > oo it is seen that, when pi; stays finite, the power transfer 
becomes independent of pj and is proportional to the number of ions in 
the upper state. The power transfer then found could be called spontaneous 
emission of phonons analogous to the situation with photons. 

For the case hay <kTz, (1.76) can be rewritten with (1.52) as: 


a 1 l 1 1 
ew Pug = ZauaTo(q, — ze) = eT 2g, — ze) 
where the power transfer constant 7a for the direct process equals 
(1.78) na = > Nita = > (heoig)*|<t|H1'(S)|9> 0? 
ij 3] 


The summation in (1.78) may be performed by a diagonal sum. 
For the case of Cr++* in CrK alum we found (ef. [1.5.4]) that « equals 
3 for the adiabatic terms in the direct process. For nq we then find: 


(1.79) na © Tr[[41' 4, |4 5, |. 
The non-adiabatic terms would require the calculation of a sixth moment. 


[1.7.2] Indirect processes. To calculate the power transfer due 
to the indirect processes, the same procedure as above is followed. For the 
transition rates we take the integrant of (1.55) again without the shape 
factor g(m@—ai). The power transferred is then found to be 


- NAN h 27] 1 
\ Pose =a iB 0(@x) e(@1)|<t|H2'|7>|? os(Hi) ae _ 7) dor,dw, = 
uj 
0 


k ah aE 
(1.80) BES Sh 
| = F(t 7) 
The power transfer constant 7,,, due to the indirect processes becomes 
(1-31) A > Nina = Y A (hovi)?| i] H2"|7>[? TF 
ij 


where A is a constant and « depends on whether adiabatic or non-adiabatic 
terms are important, and on the integration over the densities of the 
lattice vibrations (cf. (1.56)). The summation can again be performed 
by a diagonal sum 


(1.82) feoale (Hie toe | 


[1.7.3] The influence of exchange interactions on the power 
transfer for the case of S=1/2 is different for the direct and indirect 
processes. For S=1/2, #1 and #2 can only contain the spin operators 
linearly. In section [1.4] (cf. (1.33)) we have discussed that the 
fourth moment of the spin operators >;Sz,y,z, is influenced by the 
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exchange but because of commutation properties the second moment is 
not. Therefore yq will be influenced by exchange interactions. If w,, > OL 
the lattice vibrations up to the frequency w,, can be excited. 

The value of 7,,, does not depend on exchange interactions when all 
ions are equivalent. For S>1/2 each salt has to be considered separately. 


4 ‘) 
1. cr eX 


Lattice 


> 


Ee 


Fig. 1.6. Model for the energy balance. 


[1.7.4] The power balance equations are obtained by com- 
bining (1.73) and (1.77) or (1.80). The situation is sketched in fig. 1.6. 
The rate of change of the internal energy Us of the spin system is: 


dus — QW ] 1 
(1.83) di <<. ~ 1717 — 7): 


And the steady state solution: 
(1.84) Ts = Ty, + oa 
which with (1.73) results in: 


oe) Ps = Faw 


and for the case where (1.26) is valid 


1.86 ee ee 
(1.86) "Ist SnD,  2k(2S-+1)q Ts 


To solve for the transient solution of (1.83) we realize that for H kept 
constant 


(1.87) dUs _ b+ OH? dT's 
dt Ts? dt 


where (b-++-CH?)/7T's2 is the spin specific heat at constant magnetic field 
(cf. [11]), and C is the Curie constant. The term b is made up of the various 
spin interactions and the zero field splitting. 
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Solving (1.83) and (1.87) for the case that after t=0 the induced 
transition rate has the value Wand at t=0 the spin temperature is 7'5(0): 


The magnetization M=CH/7's will thus show an exponential change 
towards equilibrium with a relaxation time 


b+CH? 


(1.89) Tir — QW4nT, 


The ratio between the relaxation times as found by the steady state 
method using formula (1.26) and by the transient method with W=0 
is with respect to (1.86) and (1.89) 
(1.90) =H 36-(6 1100 4 1)2 ee 
Tigt wo 

At high magnetic fields (H?*> b/C) when the relation o=yH holds, this 
factor becomes 1, 4, 10, 20, 35 for S=1/2, 1, 3/2, 2, 5/2 respectively. 

The background for the difference between 1,, and t1,, lies in the 
derivation of the equation (1.26) which was obtained by considering 
only the relaxation probability at the saturated transition, while one 
actually has to consider the total power transfer between spin system 
and lattice as was done in formulas (1.77) and (1.83) leading to (1.85). 


[1.7.5] For the case W=0 the values of 7,, in (1.89) due to the 
direct process are calculated by Van ViEcK [36] for CrK alum and are: 
T14,=11; 9; 6.7; and 3 ms for H equal to 0; 0.5; 1 and 3 Kg respectively. 

The order of magnitude compares well with the experimental data; 
however, the temperature and field dependence for concentrated sub- 
stances are completely different. This will be discussed later. For the 
indirect processes VAN VLECK finds 


6+CH? 
(197) Tltr = 0 5 CH? 


This formula is identical with the formula used earlier by Brons [49] 
to describe the experimental results at liquid air temperatures. The value 
of p for CrK alum is expected to be between 0.6 and 0.4, dependent on 
which of the normal modes Q; of the water octahedron is the most effective 
to produce relaxation, but independent of the values assigned to the 
cubic field splitting of the orbits and the zero field spin splitting. The value 
of » depends on the commutator (1.82) and is expected also not to depend 
on temperature. For most salts this is found not to be the case, however 
(cf. [11]). This discrepancy is rather serious and must not be dismissed 
lightly. To obtain more insight into the processes involved it seems 
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advisable to determine experimentally and theoretically the anisotropy 
of p (or Mpa) and the influence of the exchange interactions. 

Van VuiEck has averaged his values of the power transfer over a 
gaussian distribution of the magnetic field, thus taking into account the 
effect. of dipolar fields. It is better, however, to take in the commutators 
for calculating 7, the total spin Hamiltonian with spin-spin interactions, 
making the calculations exact. This will be considered more fully in 
Chapter IV. 


[1.7.6] Deviations from Curie’s law. For the case that the 
spin splittings become of the order of k7's, we have to contain the density 
matrix in (1.76) and (1.80). The magnetic absorption is most easily found 
by using one of the Kramers—Kronie relations [50] [51] 


co 
: j ' 2 ¢ wx’(o) 
(1.92) x{ (on) — (00) = = [ SE do 
0 


This relation only holds when the magnetization is a linear function of 
h, so at low r.f. fields. 
For the static susceptibility one finds 


da. 


2 x(@) 
HO Noo 5 } o 
0 

When the line width of the transition is temperature independent 
we see that the temperature dependence of y” and yo are the same. 

For many salts the Curie-Weiss law holds 

ff « l > 
(1.93) x0 OC Xx T,—6 for T's 2 6. 
At higher r.f. fields one can in first approximation insert 7's instead of 
the bath temperature as is done in (1.93). 

When short range order effects are neglected, the curves yo/z" vs. h2 
as drawn in fig. 1.3 will still be straight lines, with the same slope, but 
having the ordinate (1—6/7',) at h2=0. So for an antiferromagnetic 
coupling they are higher, and for a ferromagnetic one lower. We assumed 
the line width to be temperature independent but this is not realised 
for dy we 0. 


[1.8] The relaxation method 


When the field applied to a paramagnetic substance is changed, the 
magnetization will be able to follow this change if the rate is slow enough. 
At higher rates M will have a phase lag, with subsequent absorption and 
dispersion. This provides a method for determining 1, the so-called 


relaxation method. The technique has been developed mainly by GorTER 
e.a. [11] [52]. 
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The spin system is again regarded as an ensemble in internal equilibrium 
at a temperature 7's, which was first done by Casmuir and Du Pri [53]. 
We can start therefore with (1.77) setting W=0 


(1.94) AUg= —n/T,AT dt. 


We are interested in cases where 7's ~ 7';, as the relative variations in 
# are kept small. Changing H will change the internal energy by an amount 


(1.95) ie Cu(sa7), 4M +Cm (35), AH 


where Cy and Cy are the specific heats of the spin system at constant 
field and constant magnetization respectively. 
On the other hand we have 


(1.96) AT =Ts—T, = (<5) 4M + (35), 4H. 


From equations (1.93), (1.94) and (1.95) AZ’ can be eliminated and one 
arrives at the following equation 


= ey oT AV eet) or a 
(1.97) ( ae pat) (s7),40 4 (Cm «4 dt) (Sa), 44 210) 
We now superpose on H a parallel sinusoidal component and look for 


the corresponding component in M and 7's 


{ H=H,+he 
eg, | M=M,+ mem", 
The static and the complex susceptibility are defined as resp. 


OM, £. - uw 
= (=), and ¥ = y¥'—1ty" = m/h. 


Inserting these definitions in (1.96) results in the equation for the 
frequency dependent susceptibility. 


__ ioOu+n/T 
(1.99) HO = FoGaq tent 
or with 
On -Csiamen CE? 
(1.100) Some SINISE 
nT  b+CH? 
(1.101) Lees anes D 
one finds 
F 


(1.102) VAD ha A NE aires 
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At low frequencies one finds normally %=Zo, and at high frequencies 


CG 
tea = (1—F) x0 = Ge Xo» 


the so-called adiabatic susceptibility. At these high frequencies the spin 
system is supposed to be still in internal equilibrium. 


AID Oe & 
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(To be continued ) 
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PHYSICS 


ON THE POWER TRANSFER BETWEEN PARAMAGNETIC SPINS 
AND CRYSTAL LATTICE. I 
EXPERIMENTAL METHODS AND TECHNIQUES 
Bx 


B. BOLGER 


(Communicated by Prof. C. J. GorreR at the meeting of June 27, 1959) 


[2.1] Introduction 

In this chapter we shall discuss some specific techniques for determining 
the spin lattice relaxation time by means of paramagnetic resonance. 
There are two in principle different methods, which have been used up 
till now. With the first method, the steady state saturation technique, 
the decrease in relative energy absorption is measured as the microwave 
field strength is increased. The power transfer constant 7 or the relaxation 
time t is determined by making use of (1.85) and (1.86). 

The second method displays with a low power r.f. signal the mag- 
netic absorption and its return to equilibrium after being saturated with a 
pulse of high rf. field. With the equations (1.88) or (1.65) the relaxation 
behaviour can be determined. 

Steady state saturation is a heat transfer problem and the spins have 
a forced motion. The pulse method is more analogous to the relaxation 
method; a kind of RC time is measured. During the time the high power 
pulse is off, the spins describe a free-running precession. In the pulse 
method one can choose the saturating frequency to be different from the 
sensing frequency so that one transition is saturated and its influence 
on the energy level populations of other transitions is measured. These 
methods will be discussed at the end of this chapter. 


[2.2] Magnetic losses in resonant structures 

To increase the effects of magnetic resonance the samples are placed 
in a low loss resonant structure. Maxwell’s equations for the electro- 
magnetic field inside the resonator have therefore to be combined with 
the equations for the magnetization of the sample, 

An elegant method of describing the behaviour of resonant structures 
has been given by PANNENBORG [1], by making use of the scattering 
matrix. The scattering matrix S of a junction relates the vectors of 
the incoming to those of the outgoing waves. One of its most important 
properties is that for a lossless junction S is unitary. By lumping the losses 
inside a resonator into an extra output line with a matched termination, 
its scattering matrix becomes unitary. The coupling of this line has to 
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be made variable when the magnetic losses are changed. When non-linear 
effects of the magnetization become important, such as in parametric 
amplifiers etc., one has to go back to Maxwell’s equations. 

If the modes of vibration of the resonator are well separated, one can 
make fruitful use of the equivalent circuit method, where all the distributed 
losses and stored energies are concentrated in resistances and reactances. 
One of the disadvantages of this method is that one has to assign a 
specific location to the distributed circuit parameters. At another point 
along the transmission line, transformation of the impedance Z takes 
place, due to the finite wavelength and according to the well-known 
transmission line equations [2]. Therefore all circuit representations will 
from now on refer to a certain reference plane in the transmission line. 
For the details of microwave theory we refer to the current handbooks [3]. 


Ze eeEG VR M n*Zq «jB R M 
as 
5 Pf Pi4 
a b 


Fig. 2.1. a) Equivalent circuit of resonator with magnetic material. 
b) Transformed equivalent circuit. 


The resonator without magnetic material can be represented as a 
simple series RLC circuit (cf. fig. 2.1) with a resistance R due to the losses 
in cavity walls and dielectric materials present. The reactance at frequency 
o equals 


(2.01) jB = j(2-&) QoR 


where we is the cavity resonance frequency 1//L’C and Qo=oL’'/R the 
quality factor. 

The influence of the magnetic material can be represented by the 
parallel combination of resistance M and reactance /X, due to the magnetic 
absorption and dispersion, in series with R. . 

If the paramagnetic resonance line has a Lorentzian shape we can with 
(1.13) write for the total impedance of the magnetic material, 


(2.02) Zm=M/1+7(w—ox)t2. 


The cavity is coupled to a generator and a load. This coupling can be 
represented by transformers with transformation ratios of 1:n and 
p:1 resp. The generator with voltage # at frequency is supposed to 
be matched to the input transmission line with a characteristic impedance 
Zo. The load and the generator mesh can be transformed into the cavity 


mesh (fig. 2.1b). 
Some of the power falling on the cavity will be reflected towards the 
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generator, some will be transmitted towards the load, both depending 
on Zm. This provides two ways of measuring magnetic r.f. susceptibilities. 
by transmission and by reflection. 

Quantities more accessible to physical interpretation than impedances 
are quality factors. Taking for the energy stored in the volume Ve of 
the resonator 


i 1 way Ab) 'h2 ——e 1 'h2 i 
(2.03) GS im! (e’H2 + y'h®)dV = a fh dy 
Vie Ve 
then the quality factors are defined as: 
oh. ® X energy stored - oG 
Yo = power diss. inside the cay. without magnetic abs. @R 
o oG = oG 
Qe, -_- power dissipated in the generator impedance n®? Zot? 
; aG oG 
(2.04) Qe, = power dissipated in the load impedance pli 
Q oG 7 oG 
- ~~ total power dissipated without magnetic abs. @2(n2Zo-+p2L+-R) 
Om = woh = of 
sd power dissipated in the magnetic material § 72M’ 


As the inverse of the Q’s we define u=1/Q with the same indices. As 
adding of inverse Q’s has often a simple physical meaning (adding of 
losses) we shall use mainly the w’s in the following. 

The magnetically dissipated power is according to (1.14) 


(2.05) Pm=1/2 wy" f[ h?dV im, 


integrated over the volume V» of the magnetic material and so we have 
with (2.03) and (2.04); 
(2.06) Um = G- = dof" 


where the filling factor / equals: 


(2.07) f= Jh?dV/ f uw'h2dV. 
Vin Ve 


[2.3] The reflection method 


[2.3.1] The reflection method to detect paramagnetic resonance 
and saturation will first be dealt with. The reflection coefficient I’ at a 
certain reference plane is defined as the ratio of the vector of the electric 
field of the backwards travelling wave to that of the forwards travelling 
wave, both waves in the same mode of propagation of the transmission line. 
For a transmission line with characteristic impedance Zp terminated by 
an impedance Z we have by standard microwave theory 


(2.08) pa=47% 


in the plane of Z or equally 


(2.09) 4 14? 
Zo 1-I" 


In fig. 2.1b the load impedance can be omitted and is taken zero. For 
I’ one obtains: 


(2.10) jee jB+R+M/[1+j(@ —az)t2] —n2Zo 
jB+R+M/[1+7(@—o1)t2] +n? Zo" 


Multiplying numerator and denominator by i2/oG, and using the defi- 
nitions of the w’s: 


- (@ —We) 
72 
LS ee 


Gun ca 


+o +Um{1 —j(@ —az)t2} —Ue, 


TP Uo +Um {1 —j(@ —az)T2} +-Ue, 


@e and wy are the cavity resonance and Larmor frequencies respectively. 

To determine wm two measurements are necessary. In the first place 
a determination of the reflection coefficient [%m with the magnetic field 
adjusted to the centre of the magnetic transition to be studied. Secondly 
the reflection coefficient [> with the magnetic field at a high value, making 
Um=90. This provides a determination of the relative external loading. 
When for both measurements the cavity is adjusted to resonance, then 
I and I’m are real numbers. In the following this adjustment is always 
supposed to be made. We keep the sign for I’, however, and so I is 
negative when u,>wo corresponding to an overcoupled cavity, and I is 
positive when ue<wuo for undercoupling. When wm is positive we have 
Im>TIo for both cases. 

From the above-mentioned two measurements wp, can be calculated 
as we have with (2.11); 


ae) ee eee es ey lee (lat) 


where P,,, is the power dissipated in the external loading of the generator 
impedance, while there is also magnetic absorption. 

Of the power Pi, incident on the iris, the part P,,, dissipated inside the 
cavity is, whilst there is also magnetic absorption, 


Ve 1L— I'm?) Pi; 


cay = | 
and the reflected power is 

(2.13) Pra Ly Ps. 

The power dissipated externally can be found from (2.11) if o=@c=or1 


to be: 
(2.14) (en See fad i I 


Uo +Um 
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If HE is the electric field strength of the incident wave, then (1 —I\E is 
the electric field inside the cavity. As the external loading does not change 
during the two experiments discussed above, we have with (2.14), (2.03) 
and (2.04) that P,,,« W&«h2. The proportionality factor between P.x¢ and 
h2 is only a geometrical one. 
When the rf. field is uniform over the sample dimensions the final 


form of the equations becomes: 


et Tr, Pext | Pext 


(2.15) Pn QWy 17 
¢ a Pext _ (1—I'm) (1 —L0) 
(2.16) (ghee 3) PST 
(2.17) Pj =(l—Palhs 


[2.3.2] Determination of 7. A fast method of determining the 
reflection coefficients is to measure the power P, reflected by the cavity 
with a microwave bridge. If P; is the incident power, J” follows from (2.13). 

The first term of the right hand side of (2.15) is independent of h 
because Wi; « P,,, * 42, giving a straight line for the P..,/Pm & 1/z" versus 
Px plot if 7 and 7’, are independent of h. The slope is 1/7. 

To calculate t with (1.86) or (1.101) one only has to know, moreover, 
the weight of the sample and the frequency o. 

For the case of inhomogeneously broadened resonance lines the 
saturation behaviour becomes different (cf. [{1.4.4]). 

The determination of 1 by the method of BLOEMBERGEN, PURCELL and 
Pounp reduces to measuring the r.f. field hy, at which 7” has dropped to 
half its initial value. This occurs when the first and second terms on the 
right side of (2.15) become equal so that »=QWy/Tx or with (1.26) 


: ; . 1 
(2.18) 1/2 y®A¥, 9(0)|<i|Salj>|2 = 2. 
This method involves a calibration of h, — which is difficult to obtain 
more precise than 30%, — and knowledge of g(0) the matrix element. 


A serious disadvantage of Bloembergen’s method is that when 7 or Tz 
depend on h, of which we shall see examples later, one does not use all 
the information available in the measurements. 


[2.3.3] The polar impedance diagram. We shall briefly sketch 
what happens in the polar impedance diagram (cf. [2]). This diagram or 
Smith chart (fig. 2.2) represents the complex I’ plane. The loci of constant 
resistance and constant reactance form two families of orthogonal circles. 
The constant reactance circles, dotted-dashed in fig. 2.2, have centres on 
Re I7=1; the constant resistance ones, the solid lines in fig. 2.2, their 
centres on the real I" axis. 

When varying the frequency of the generator, the impedance presented 
by an empty resonator without any other resonances nearby will, at 
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a certain reference plane, be a locus of constant resistance, as only the 
Teactance is varied. The reflection coefficient is the vector drawn from 
the origin to the locus. Curve (a) corresponds to total reflection, R—=0 
curve (b) R<Zp corresponds to an overcoupled cavity, curve (c) Raa 
to a matched cavity and curve (d) R>Zp to an undercoupled cavity. 


Fig. 2.2. Polar impedance diagram. 


resistance circles. ape 0) 
eae reactance circles. bei Fo 
——-— magnetic effect, H varied. c R=Zp 

ad Zo 


Suppose now the cavity is at resonance (J/’ real), then varying wy 
(H) will change Z according to the dashed circles due to the influence of 
the magnetic material. The reactance part is due to dispersion, the 
resistance to absorption. One observes that magnetic absorption increases 
the reflected power for undercoupled and matched cavities and decreases 
it for overcoupled cavities. As long as the magnetic material is dissipative, 
its impedance locus has of course to stay inside the cavity locus. Magnetic 
materials can be brought into an emissive state. Their locus is then 
inverted with the cavity resonance point as the centre of inversion. 


[2.4] The transmission method 


The calculating labour involved to obtain the saturation curves with 
the reflection method is quite substantial and we therefore tried to find 
another method which required less calculating time. We should like a 
formula like (2.15) where instead of P,,, a variable is used, proportional 
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to the square of the radio frequency field and directly accessible to 
measurement. Inspecting fig. 2.1 one sees that by inserting for the load 
a square law detector with detector sensitivity D the signal « becomes: 


(2.19) Du=P,=iL he. 


If zm is the signal whilst the cavity and magnetization are in resonance, 
one finds with (2.04) and (2.11): 


4te, Ue, P; 


(2.20) va om = Dig = (wr +-tm)?2 


and without magnetic absorption for the signal Xo: 


4Ue, U 
(2.21) Pr = Da Ss oe 


Dividing (2.21) by (2.20) gives 


For the absorbed power we can now write with (2.04) 
= Oe Dae, ws Mt wo 
(2.22) Pm = Ditm = & Dim (2 1), 


In these formulae there are two unknown external couplings, one can 
be eliminated, however, by varying the generator coupling with fixed 
load coupling so as to maximize Py (2.21). This makes up=2ue, and 


(2.23) Py = % D = Py, 


Together with (2.22) follows 


ep spe 4B)" 


which is the formula to be used for calculating the saturation graphs. 
Pi/v is an easily calibrated constant and am«h2. If the difference 
xX —Xm=dKxy we may expand the root: 


(2 25) Xm P;/ xo vm 
; Pu le 


As 6 and xm are directly accessible to measurement, the calculations take 


very little time. The power dissipated by the magnetic material is then 
found to be 


(2.26) Foe ces. 
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[2.5] Apparatus for the reflection measurements 


The experiments on the steady state saturation have been carried 
out at microwave frequencies around 10 KMHz. For the general theory 
and techniques of microwaves reference is made to the current handbooks 
on the subject, especially those of the M.I.T. series nos. 8, 9, 10 and 11. 
Most of the microwave plumbing such as attenuators, tuners, cavities ete. 
have been designed by us and were made in one of the laboratory work- 
shops. By measuring the power incident on and reflected by a cavity, the 
reflection coefficient can be obtained from P,=I72P;. To measure these 
two powers separately, use is made of a microwave bridge (fig. 2.3). The 
bridge element was a 3db four outlet directional coupler of the RiBLer 
and Saap type [4]. This element consists of pairs of coupling slots and 
is by far superior to the usual magic Tee because of the inherent directivity 
of each pair of slots. 


Fig. 2.3. Apparatus for measurements of reflection coefficients. 


B,; and Bz barrettors; Vi... V«6 attenuators; L matched load; M vacuum seal; 
S' sample in 2/4, long cavity. 


Initially a magic Tee has been used but sliding screw tuners were 
needed to obtain satisfactory behaviour of this bridge element. 

For the high power microwave source was used a Philips 55395, a two 
cavity klystron, able to deliver 140 to 200 Watt C.W. power in the 
X band [5]. 

Because of the high powers used, special care has to be taken to prevent 
burning out of attenuators and detectors. For V3 we constructed a water 
attenuator, while the fixed attenuators Va, Vs and Vg were made of a 
ceramic caslode. Their size was chosen so as to keep the power dissipation in 
the calibrated variable attenuators Vi and V2 below the burning out level. 
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These variable attenuators contained a platinized glass vane as dissipative 
element. They were calibrated 3 times at three frequencies far apart with 
the aid of a barrettor bridge and also compared to each other. The repro- 
ducibility was within 0.1 db. 

The matched terminations were made from double tapered caslode 
and had a voltage standing wave ratio below 1.1 in the frequency range 
used. A handy sliding match was made from a long piece of tapered oak 
just fitting inside the wave guide (V.S.W.R.< 1.08). 

For power detection we used silicon crystals or barrettors. Barrettors 
contain as sensitive element a very thin platinum wire which changes 
its resistance upon heating. Power measurements are made by using the 
equivalence of heating by microwaves or by a direct current. The barrettors 
were of the 821 type (Sperry or Narda Co.) with an operating resistance 
of 2002. For power detection, barrettors compare favourably with 
crystals because of their good square law and burning out characteristics. 
They were inserted in one branch of a Wheatstone bridge, which could 
be balanced by varying the bridge current, causing a change in the 
barrettor resistance B. This current was measured by a Tinsley potentio- 
meter. The microwave power can be found from 


(2.27) P= Blig?—i2). 


i and ip are the balancing currents required with and without microwave 
power respectively. 

As we are interested in a change in the cavity reflection or transmission, 
due to saturation of the magnetic losses, when varying the incident 
power, a null method was used. Two barrettors were placed in neigh- 
bouring branches of a Wheatstone bridge (fig. 2.4), passing the same 
d.c. current. In all measurements the bridge is balanced to zero galvano- 
meter reading and the bridge current measured by the voltage drop over 
the 1 2 standard resistance by a Tinsley potentiometer. As the barrettors 
may need slightly different heating powers one could either give one 
barrettor an additional heating, by for instance a 100 KHz signal, or 
take the following procedure. Set switches S; and S: to the 200 Q standard 
resistances A; and Rs, adjust G to zero by R5. Set S; to barrettor By, and 
balance with bridge current by R7 or Rs. This makes By; equal to 200 Q. 
Set S: to By and balance with Rs keeping the bridge current the same. 
The bridge is now ready for power measurements. The incident microwave 
power P; is measured by J, (in fig. 2.3), when setting S; to B,; and Sp» 
to Ry and balancing by changing the bridge current. The difference between 
P, and P; is measured by setting S; to By and S, to Bs and balancing with 
the attenuator V»2 in fig. 2.3. 

A calibration for the difference in attenuation in the microwave bridge 
arms and the difference in the detecting sensitivities of B, and Bs is 
obtained by replacing the cavity with a total reflection. To take measure- 
ments in quick succession it is necessary to heat the barrettors also to 
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about 200 2 when not passing the bridge current. The time constant 
involved is of the order of one or two minutes and not at all the generally 
used time constant characterizing the ability of a barrettor to follow a 
modulation of the microwave power, which is of the order of 1 ms. The 
time constant involved here is much longer as it is determined by the 
establishment of a temperature distribution along the barrettor wire and 
holder after changing the wire dissipation. 


[ae 
Ry Ra. ya) 
LALA =. y: 


Zi TINSLEY Pot. 


Fig. 2.4. Barrettor bridge. B, and By barrettors. 


It is more difficult to adjust the frequency of a two-cavity klystron 
than of a reflex klystron. Therefore the cavities had to be made tunable. 
Two different cavities with their cryostats were used. For the first cryostat 
a german silver waveguide, silvered inside, went down into the liquid 
helium bath, vacuum sealed on top by two thin mica windows 1/4 dg apart. 
At the end was a rectangular cavity 24, long with an adjustable iris. 
The cavity was tunable by a plunger movable from the top of the cryostat. 
The coupling could be adjusted before cooling to the desired value by 
placing a small capacitive screw in the middle of the iris. By screwing 
it in one initially increases the coupling, which decreases again after the 
iris has become resonant. 

As it is not possible to adjust this coupling during a run and as some- 
times the plunger stuck because of freezing, a different cavity of about 
80 cm, or about 18 A, long was constructed. This cavity was tunable 
with a teflon flap and the iris consisted of a movable copper sheet. Iris 
and tuner were located above the vacuum seal. Since the filling factor 
of the paramagnetic sample decreased more than the @ increased, there 
was a loss in sensitivity, but having ample sensitivity to start with we 
could spare some in exchange for convenience. The movable iris in this 
case provided, when closed, the total reflection needed for calibration 
of the barrettor sensitivities as discussed before. For the first cavity [’= 1 
was obtained by detuning the cavity far from resonance. That this really 
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provided a total reflection had been checked before by interchanging the 
detuned cavity by a total reflection. The unloaded cavity Q’s were of 
the order of 2800. This low value of Qo was taken so as to make the whole 
set up not too frequency-sensitive. To prevent boiling bubbles of the 
helium inside the cavity to change the resonant frequency, this one 
was filled with foamplastic. The low temperature techniques used were 
those customary in our laboratory. The magnet had a resistance of 2 a 
and provided a field of 10 Ko at 100 A current within a gap of about 
9em and was constructed in the laboratory [6]. 


[2.6] Apparatus for the transmission measurements 


The experiments with the transmission method were also carried out 
at microwave frequencies around 10 KMHz. The apparatus is sketched 
in fig. 2.5. 


V. U 
kiystron = 


Fig. 2.5. Apparatus for transmission measurements. 


B directional coupler; Ji, J, input, output coupling; ZL matched load; My, Me 
vacuum seals; S sample; U unilines; V attenuators; Xj, X2 detecting crystals. 


The same power source was used as with the reflection measurements. 
The power was detected by two carefully selected silicon crystals. 
The requirements for the selection were: a very close to square law 
detection and a difference in output between the two of less than 0.1 mV 
in the range from 0 to 20 mV, the difference being made zero at 20 mV. 

In fig. 2.5 the incident power is sampled by crystal X1, through a 3 db 
directional coupler [2.5] and the transmitted power by Xo. Their output 
voltages will be called a and a: respectively from now on. The outputs 
were always kept below or equal to 20 mV and were measured by a 
Philips d.e. millivoltmeter G.M. 6010. 
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The cavity was about 80cm long and tunable by a dielectric vane of 
teflon (Qo ~ 2800). The irises had to be made variable, the input iris J, 
in order to match the cavity to the generator, a condition set to the 
validity of equation (2.24), the output iris Jz in order to be able to work 
in a convenient range of V2. 

An easy way of doing this is by a movable shorting plunger in the 
output waveguide [7] as shown in fig. 2.5. The coupling takes place 
through the magnetic h.f. fields and for iris 7; for instance, maximum 
coupling is obtained when the short is 1/4 dg from the centre of the iris, 
minimum at 1/. A, distance. For J. the situation is reversed. The two 
irises were placed before the tuner so that their electrical distance did 
not change when the insertion of the dielectric was altered. 

For a series of measurements the following procedure was used. 

At a static magnetic field, chosen high enough to make the magnetic 
losses zero, the coupling of iris 7; was adjusted so as to maximize the 
transmission of the tuned cavity, Jz being chosen beforehand to a con- 
venient value. The microwave power is now set to make z2=20 mV and 
the difference 6 between 2; and a2 made zero, by changing V;. During the 
rest of this series of measurements the attenuator setting is not changed. 

The transmitted power without magnetic absorption is given by 41, 
with magnetic absorption by x2. They correspond to xp and xm respectivily 
in formulae (2.24) and (2.25) and can be simultaneously measured. 
This can most conveniently be done by a X—Y recorder. Most of the 
time we measured zv2—x%1=6 and 2, as they enter directly into equation 
(2.25). When the tracking of X; and X2 was not according to the require- 
ments set, a correction was made. The value of the detector sensitivity 
P;/x; in mW/mV is calibrated by substituting a barrettor for the cavity 
at the input branch. If this sensitivity is C db, then the input power 
during a run is P;=C+V,+10 log a dbm, db above the 1 mW level. 
Vi is the attenuator setting in db. 


[2.7] Pulse methods for determining the relaxation times 


The principle of this method has been indicated in [2.1]. The approach 
of the magnetization towards its steady state value during or after a 
pulse of saturating power will depend on the transition observed and 
the initial situation. For the case of non-important cross relaxations, the 
considerations in [1.6.2] provide the necessary relations for interpreting 
the observations. For the case of internal spin equilibrium (cf. (1.88)), or 
one of the time constants 1/A, (cf. (1.65)) much larger than the others, the 
change of the magnetization is according to a simple exponential. 

The magnetization can be observed by a microwave signal of frequency 
ws near the transition to be studied and an intensity so low as not to 
disturb the system appreciably. We will call this microwave field the 
sensing or signal field. The magnetization can be made to deviate from 
its equilibrium value by a pulse of high microwave power, called the 
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pump power, at frequency wy. The pump may saturate the signal tran- 
sition or one of the others when present. Observation of the fall of the 
magnetization during the high power pulse may provide information 
about the cross relaxations. 

The situation where wp—«s has been employed by BowErs and Mims [8] 
and also by Davis and StRanDBERG [9]. An apparatus for this type of 
measurements is under test in our group, but only preliminary measure- 
ments have been made with it. In our set up @s/272 can be changed in 
steps of about 100 MHz to other cavity resonances in the range of 8.3 
to 11 KMHz. The apparatus is essentially the same as used for the trans- 
mission measurements. The signal frequency is generated by a low power 
klystron and is detected by a superheterodyne receiver with a sensitivity 
of about 10-12 W. 

Experiments with ws/227 = 1420 MHz and w/2a= 4000 MHz or 10 KMHz 
were made in connection with the construction of a low noise microwave 
amplifier. From the characteristics of some types of these amplifiers it 
is possible to obtain information about the relaxation processes. We shall 
therefore mention briefly the principle of a so called three level maser 
as constructed by our group. 


[2.8] The principle of operation of a three level maser 


[2.8.1] Introduction. In conventional microwave amplifiers (travel- 
ling wave tubes and klystron amplifiers) the source of energy is the 
kinetic energy of an electron beam which introduces a certain amount 
of noise into the amplifier. Recently it has been found possible to amplify 
a microwave signal by converting the internal energy of an atomic system 
directly into microwave energy. When the equilibrium of the system is 
disturbed in a particular manner it is possible to obtain an emissive 
condition, and when the medium is stimulated by a microwave signal, 
the energy thus stored can be emitted coherently. GoRDON, ZEIGER and 
TowNnEs [10] have coined the term “‘maser’’ to describe such “microwave 
amplification by stimulated emission of radiation”. The small probability 
of spontaneous emission at radio frequencies (the main source of noise) 
makes it possible for such an amplifier to have a very low noise figure. 

Of the many types of masers (discussed by WrrrKe [11]) the only one 
giving continuous emission is the three level maser (BLOEMBERGEN [12]) 
using a paramagnetic material as the working substance. 

For a detailed analysis of the characteristics of three level masers 
reference is made to the papers by BurcHEr [13] and for the work on 
this subject by our group to [14]. 


[2.8.2] Principle of operation. In section [1.6] we discussed the 
influence of microwave fields on the energy level populations of a multi- 
level spin system with negligible cross relaxations. Inspecting formula 
(1.70) it is seen that the difference in population Aj, of level j and k can 
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become negative when Cra,12 > Cyi,i2 and Wiz (at the pump frequency) 
is large enough, while levels 1 and 2 have to be non-adjacent ones. 
The power absorbed by the j to k transition also becomes negative 
when Aj, becomes negative (cf. (1.71)) or in other words we have attained 
emission. 

The frequency of the 7 to & transition will be taken to be the signal 
frequency ws. Suppose the sample is placed in a reflection cavity. When 
the frequencies in formula (2.11) are so adjusted that WO=WL=s, it is 
seen that the reflected power becomes larger than the incident power 
when —wUm> wo (I’m<—1) so when the emitted power is larger than that 
dissipated in the cavity walls (w=1/Q were defined in (2.04)). For this 
case we have obtained amplification. The bandwidth can be obtained 
from the total quality factor (Qm is negative). Under certain conditions 
the power emitted by the sample will be so large that the amplifier breaks 
into oscillation. 

Under the conditions of oscillation the signal intensity will build up, 
until the transition is so far saturated that the power emitted by the 
sample just compensates for the losses in the cavity and external coupling. 
The level of oscillation can be obtained from this condition 


(2.28) Um + Uo + Ue = 0. 


In (1.70) Ajx is determined by the pump intensity for small signal 

powers. For large signals, however, 4; will also depend on its intensity 
due to the presence of Vj,~=Ujx4-+ Wjx in D and C’s in (1.69). 
For a spin system with only three levels calculation of (1.70) becomes 
easy. Numerating the levels from the lowest and pumping between levels 
1 and 3 results for high pump powers (W143 large) in an induced absorption 
at M23 of 


32 U 23 — 21 U 12 


Ue1+U 23+ Wes 


Nh?o32 
(2.29) Pog = Sr Wes 


where 2Uy7= Uyt U jt. 
The condition for obtaining negative values of P23 is seen to be 


(2.30) 21042 > 032 93. 


At low signal powers (W23 < Us: + U2s) the ratio of the emitted power 
when Wy3 co to the absorbed power when Wi3=0 is, when Wog is 
the same for both cases: 


Pen Umem 21 Ui2—32U 28 


(2.31) Pra,  Umap w32( Ue1 + Ues) 


which depends only on the ratios of relaxation rates and frequencies of 
transitions 23 and 12, but not on filling factors, cavity Q’s, sample weight 
etc. In case of oscillations, we have mostly We3 >> Ue + Uo3. When this 
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condition is fulfilled over the dimension of the sample the maximum 


power emitted becomes, still with Wi3 > oo; 
N h?ag32 aa Ty 
(2.32) P28max = 3 AT (21 U2 — 32 U3). 
It is noted that this expression does neither depend on line width, 
when homogeneously broadened, nor on matrix element and filling factor. 


[2.8.3] In the experimental arrangement the crystals were 
placed in cavities resonant simultaneously at the signal and pump 
frequencies, and cooled by a refrigerant in a Dewar-glass. For the best 
amplifying characteristics the Qo of the cavity must be high at the signal 
frequency and preferably also at the pump frequency. The filling factor 
} for the signal frequency must be as high as is consistent with the require- 
ment to saturate throughout the crystal with the pump field. 

A number of cavities has been used to obtain a variety of field con- 
figurations for the different frequencies, needed for the substances tested. 
Cavities used for 1420 and 3850 MHz were: 

(1) A coaxial cavity 2/2 long at 1420 MHz with a 7'#y11 mode at 
3850 MHz, and the 34/2 coaxial mode for 4260 MHz. This cavity had a 
high Qo but the filling factor was poor. To increase / we tried: 

(2) A double reentrant cavity approximately 4/4 at 1420 MHz and 3/4/ 
at 3850 MHz. To obtain the right electrical length for both frequencies, a 
dise of teflon or quartz was inserted at the electrical antinode of the 
3850 MHz mode. As the 1420 MHz electric field was small in that region 
the disc mainly influenced the pump field pattern. Qo was 4000 at 1.3 °K. 
Some unwanted modes were difficult to eliminate, and also because of 
another experiment to be described, the following cavity was also used. 

(3) A resonant strip 4/2 long at 1420 MHz placed in a rectangular 
cavity with either a 7’Ey03 or T’Eyo4 resonance near 3850 MHz. The strip 
was supported by two rectangular quartz rods which also shortened the 
Ag of the 3850 MHz mode. 

The cavity used for 1420 MHz and the pump frequency between 8500 
and 10700 MHz was: 

(4) A resonant strip 4/4 long at 1420 MHz at the base of a long 
rectangular cavity as described in [2.6] (ef. fig. 2.5), carrying a T’Ejo, mode. 


In each of the cavities the Qo at 1420 MHz was between 5000 and 
7000 at room temperature and the filling factors were about 30%, except 
for the first one. 


[2.8.4] The reflection characteristics were measured with a 
microwave bridge. The arrangement is sketched in fig. 2.6. The determi- 
nation of cavity quality factors and amplifier characteristics is most 
easily performed by using a frequency swept source. This was not 
available in the laboratory, so we mixed the output of a sweep generator 
(f~50 MHz) with that of a signal generator (f= 1370 MHz). Ample 
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signal strength was obtained in this way. The frequency calibration was 
obtained from markers provided by the sweep generator. For the bridge 
element a hybrid ring was constructed as described in [15]. The two 
cross branches with a characteristic impedance of (2 Zo have to be 
slightly shorter than A,/4 (because of fringing fields). The isolation 
obtained was 30 db at 1420 MHz. Cavity number 3 is sketched in figure 2.6. 
Detection took place by a superheterodyne receiver (noise figure 9 db) 
preceded by a low pass filter F (cut-off 2000 MHz); 


SUPER Ew 


LY 
4000 MHz 
OSC, 


Fig. 2.6. Apparatus for experiments with two different frequencies. B_ bridge 
element for 1420 MHz; F filter; T’ Magic Tee for 4000 MHz. 


The 4000 MHz source was a triode oscillator (EC 57) which could be 
pulse modulated. The rest of the apparatus at this frequency was a 
straightforward magic Tee bridge. 

The 10 KMHz apparatus has been discussed before in [2.6]. The klystron 
power supply was slightly changed so as to allow pulse modulation. 


[2.9] The temperature behaviour of the quality factors 


The temperature behaviour of the cavity quality factors was peculiar. 
Fig. 2.7 shows the unloaded Qo of a reentrant test cavity, measured as 
a function of temperature. This cavity was silver plated on the inside. 
The same behaviour was observed on normally polished copper cavities. 
Decreasing the temperature from 300 K on, uo behaves as expected and 
decreases also. Below the temperature where the electron mean free path. 
becomes larger than the skin depth, Preparp [16] has shown that the 
uo has to stay constant. This temperature is about 70° K for pure annealed 
copper at 1420 MHz. The increase of uo at lower temperatures to about 
the same value as at 300° K is most unexpected. We should like to explain 
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this behaviour by supposing that the surface irregularities become larger 
than the skin depth at lower temperatures, thus increasing the losses. 
A check on this hypothesis was made by polishing a copper cavity with 
the finest emery available (0.4). The increase of wo was not present 
anymore and it stayed constant below 20° K (dashed line in fig. 2.7). 

Of the cavities tried at low temperatures those made of oxygen free 
electrolytic copper (no pores) with highly polished surfaces, and with 
silver soldered joints, had the highest quality factors. 


| 
—— 


it. 4 to, © 4G) oO Gaooe 
Fig. 2.7. Qo as a function of temperature. 


[2.10] Crystal growing 


The small crystals were obtained by normal evaporation in glass discs. 
In preparing diluted crystals one has to pay attention to the fact that 
the dilution in the solid state will usually be different from that of the 
solution from which the crystal was grown. For CrK alum for instance 
the ratio of the two dilutions may be a factor 4. 

For the larger crystals, like those used for the masers, a more sophisti- 
cated growing procedure has to be followed. The crystals in question 
were Ks(Co(CN)g with a small content of KsCr(CN)¢g (~ 1 %). This last 
substance is rather unstable and the solution may not be heated too much. 

The crystals, hanging in the solutions, were rotated back and forth 
while cooling slowly in a thermostat from 32°C to about 25° C in one 
week. Rotating the crystals in one direction only, decreased the growing 
speed considerably. The solutions (without the crystals) had been saturated 
by placing them in a separate thermostat at about 35° C for some time. 
The difference between the starting temperature (32° C) and the solution 
preparation temperature should not be too large, the oversaturation will 
otherwise become too high and the solution will crystallize out immedi- 


ately. Crystals weighing about 100 g. were obtained in this way without 
too much difficulty. 


[2.11] rrors 


The attenuators used were calibrated to within 0.1 db at the frequencies 
used. The detecting sensitivities of the crystals were calibrated by using 
a barrettor, the voltage standing wave ratio of which was measured. From 
repeated measurements the accuracy of this calibration is estimated to 
be of the order of 0.1 db. 


When using the reflection method the value of J > was determined for 
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every point separately by applying a high magnetic field (10 Ko). 
The J chosen for calculating (2.16) was the average of these values. 
To prevent boiling of the liquid helium to upset the measurement, the 
cavity and waveguides inside the Dewar were filled with foam plastic. 
Nevertheless the changing helium level impaired the accuracy. The 
resulting accuracy of 7 or 1, from these measurements is estimated to 
be within 30 %. 

In the transmission method more effort was put into eliminating the 
errors, and the reproducebility of the later measurements (cf. CoNH, 
tutton salt) was much better. The main sources of errors were the drift 
of the cavity resonance, partly eliminated by taking a low Qo, the 
drift of the magnetic field and again the changing level of the liquid helium. 

The difference in detecting law of X; and X» was corrected for. The 
accuracy of the later measurements is estimated to be about 15 %. 

The pulse measurements had to be made with low signal powers and 
thus noise appeared on the pictures. If one takes the upper edge of the 
observed trace the noise does not necessarily impair the accuracy of the 
measurements. It is not excluded, however, that not all factors determining 
the time behaviour of the magnetization (such as pulse length, pump 
power etc.) have been controlled sufficiently and so some of the results 
have only a provisional character. 

(To be continued) 


REFERENCES 


PANNENBORG, A. F., Thesis, Delft (1952). 
Bronwe tt, A. B. and R. E. Beam, Theory and Application of Microwaves. 
McGraw-Hill, New York, (1947) ch. 8. 

3. M.I.T. Radiation Lab. series. Vol. 8, 9, 10 and 11. McGraw—Hill, New York. 

4, Risitet, H. J. and T. S. Saap, Proc. I.R.E. 36, 61 (1948). 

5. I am very much indebted to the N.Y. Pures Gloeilampenfabrieken for the 
generous loan of this klystron and its power supply. 

6. This magnet was designed by Prof. Dr K. W. Taconts. 

7. I am very indebted to Ir H. G. Betsers for suggesting this solution. 

8. Mins, W. B. and K. D. Bowrrs, Kam. Onnes Conf. 1958. Physica 24, 166 (1958). 

9. Davis, C. F. and M. W. P. StranpBERG, Bull. Am. Phys. Soc. Ser. II 2, 226 
(1957). 

10. Gorpon, J. P., H. J. Zeiger and C. H. Towness, Phys. Rev. 99, 1264 (1955). 

ll. Wrrrkcn, J. P., Proc. I.R.E. 45, 291 (1957). 

12. BropmMBEeRGEN, N., Phys. Rev. 104, 324 (1956). 

13. Butoner, P. N., Proc. I.R.E. 105 B (1958). 

14. Boucer, B., B. J. Rospryson and J. Ussryx, to be published in Physica. 

15. Reep, J. and G. J. Wueeter, Proc. I.R.E. Trans M.T.T. 4, 246 (1956). 

16. Preparp, A. B., Proc. Roy. Soc. A191, 385 (1947). 


Noe 


CONTENTS 


Chemistry, Physical 


Borr, J. H. pg and R. J. A. M. vAN DER Bora: The kinetics of the selective 
formation of the intermediate product in two consecutive reactions, p. 308. 

Bora, R. J. A. M. vAN DER: A comparison of an empirical and of a kinetic 
equation, both describing a selective catalytic process. (Communicated by 
Prof. J. H. pz Borr), p. 299. 

Orrens, E., A. J. Exanp and W. G. Burcers: Influence of the irradiation 
temperature on the position of the F-band in mix-crystals of KCl and 
NaCl, p. 268. 

Ortrns, E., A. J. Evanp, H. B. ZeepisK and W. G. BurcEers: Absorption 
band at 420 my in X-ray coloured KCl-crystals and in mix-crystals of 
KCl and NaCl, p. 277. 


Crystallography 
MacGILLAvRY, CAROLINE H.: Crystallographic studies on some “‘Coenzymes Q”’, 
p- 263. 


Geology 


Suat, A. F. J. and J. E. J. M. van LanpEewisE: An occurrence of corundum 
in the Precambrian of Ghana. (Communicated by Prof. G. H. R. von 
KOENIGSWALD), p. 283. 


Mechanics 


SPARENBERG, J. A.: Application of lifting surface theory to ship screws. 
(Communicated by Prof. W. P. A. vAN LAMMEREN), p. 286. 


Physics 
Bétcrr, B.: On the power transfer between paramagnetic spins and crystal 
lattice. Ia. (Communicated by Prof. C. J. Gorter), p. 315. 
Bétcrr, B.: On the power transfer between paramagnetic spins and crystal 
lattice. In. (Communicated by Prof. C. J. Gorter), p. 329. 
Béxcer, B.: On the power transfer between paramagnetic spins and crystal 
lattice. II. (Communicated by Prof. C. J. GortER), p. 348. 


PRINTED IN THE NETHERLANDS 
DRUEKERIJ HOLLAND N.V., AMSTERDAM 


